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Fig.  1.1.  Geometry  of  two  Cartesian  cross  sections,  namely  xz  and  yz,  for  the  GB 
z  =  f(x,y)  with  GBI  close  to  the  xy  plane  . 

Fig.  1.2.  GB  cross-sections:  (a)  cross  section  showing  MUs  neighboring  the  M0, 
the  best  linear  fit  with  the  normal  direction  projection;  and  (b)  cross 
section  including  single  MU. 

Fig.  1.3.  Cross-sections  generated  using  test  model  which  represents  accuracy  of 
the  GBI  estimation:  (a)  exact  result  for  R=50MU;  (b),  (c)  error  origin 
for  large  and  small  GB  curvature  radius,  respectively. 

Fig.  1.4.  GB  curvature  radius  for  most  microstructures  is  larger  than  average 
grain  size. 

Fig.  1.5.  GB  effective  mobility  v.s.  GBI  [3]:  (a)  Optical  microscope  image  of  GB; 

(b)  Dependence  of  the  relative  GB  energy  on  crystallographic 
parameters  of  the  GB  plane  in  two  adjacent  grains  (I  and  II).  The  letters 
A-FI  refer  to  eight  positions  along  GB  and  (c)  Dependence  of  the 
relative  GB  energy  on  the  tilt/twist  component  (TTC)  and 
disorientation  angle;  d)  averaged  effective  GB  mobility  as  TTC 
function  for  30-45  deg  misorientation  interval. 

Fig.  1 .6.  Cross-sections  generated  using  test  model  which  represent  initially 
spherical  grain  shrinkage  due  to  relative  GB  energy  dependence  on 
GBI  (TTC):  (a)  initial  state  R=100MU;  (b),  (c)  grain  shape  after  50  and 
100MCS.  Deviation  from  the  spherical  shape  is  evident. 

Fig.  1.7.  Modeling  volume  cross-sections:  (a)  cross  section  showing  GB  normal 
estimated  at  M0  (projection  within  the  section  plane;  and  (b)  the  same 
cross-sections  with  several  more  GB  normal  lines  projections,  at  which 
areas  excluded  from  GB  normal  estimation  are  shown  by  deep  blue 
color. 

Fig.  1.8.  (a)  Grain  shape  approximated  with  the  Dodecahedron  and  (b)  estimation 
of  GB  volume  part  for  which  GBI  can  be  derived. 

Fig.  1.9.  Cross-sections  generated  using  model  of  normal  grain  growth  with  only 
GBI  dependence  of  GB  effective  energy:  (a)  initial  microstructure;  (b, 
c)  after  20  and  200  MCS  correspondingly. 

Fig.  1.10.  Effective  GB  mobility  as  function  of  TTC  component  (a)  and  of  scalar 
GB  misorientation  (b). 

Fig.  1.11.  Cross-sections  generated  using  model  of  normal  grain  growth:  (a) 
initial  microstructure;  (b)  with  no  misorientational  but  only  GBI 
dependence  of  GB  effective  energy  after  200MCS;  and  (c,d)  with  both 
misorientational  and  GBI  dependence  of  GB  effective  energy  after  200 
and  800  MCS,  correspondingly. 

Fig.  1.12.  Effective  GB  mobility  as  function  of  scalar  GB  misorientation  only  (a); 
of  TTC  component  (b)  and  both  misorientation  and  TTC  component 

(c) . 

Fig.  1.13.  Initial  microstructure  (a)  and  initial  grain  size  distribution  (b)  for  Cases 
A  to  C. 
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Fig.  1.14.  Predicted  grain  growth  kinetics  for  Cases  A  to  C. 

Fig.  1.1 5.  Initial  microstructure  (a)  and  grain  size  distribution  (b)  for  Case  D. 

Fig.  1.16.  Schematically:  difference  in  GBI  dispersion  for  Cases  A  to  C  and  D. 

Fig.  1.17.  Predicted  grain  growth  kinetics  for  Case  D  in  comparison  with  the 
Case  C  (dashed  line). 

Fig.  1.18.  Predicted  microstructures  for  Case  D  (a)  initial,  after  150  MCS  (b)  and 
500  MCS  (c).  (d,  e).  Initial  (d)  and  predicted  grain  size  distribution 
(500  MCS). 

Fig.  1.1 9.  (a)  3D  image  of  selected  deformed  grains  and  some  rex.  nuclei,  and  (b) 
cross-section  (b)  generated  for  simulation  initial  microstructure. 

Fig.  1.20.  Stored  energy  distribution  shown  with  color  (blue  to  red  scale)  within 
the  modeling  domain.  Cross-sections  generated  transverse  plane  (a) 
and  rolling  plane  (b)  and  3D  image  (c)  showing  some  rex.  nuclei  and 
not  recrystallized  volume  after  50  MCS  (Case  A). 

Fig.  1.21.  Effective  GB  mobility  as  function  of  both  scalar  GB  misorientation 
TTC  component. 

Fig.  1.22.  MC  predicted  microstructures  after  75  MCS:  (a)  Case  A  (90%  rex.) 
and  (b)  Case  B  (80%  rex.) 

Fig.  1.23.  Predicted  recrystallization  kinetics  for  Cases  A  and  B. 

Fig.  1.24.  Predicted  microstructures  for  Cases  A  and  B  after  recrystallization 
completes:  a)  Case  A,  100  MCS  b)  Case  B,  250  MCS. 

Fig.  2.1.  Interface  for  introducing  the  spatial  texture  inhomogeneities  into  the 
modeling  volume. 

Fig.  2.2.  Local  textures  implemented  into  the  model  of  deformed  material 

microstructure.  Spatial  region  A  corresponds  to  the  001  fiber  texture, 
region  B  is  close  to  01 1  texture  component  and  region  C  is  random 
background  (undefined). 

Fig.  2.3.  Local  textures  implemented  into  the  model  of  material  microstructure, 
with  bands  consisting  from  small  subgrains  of  low  misorientation. 

Fig.  2.4.  Cross  sections  of  locally  textured  initial  modeling  volumes:  (a)  and 
(b)  BW/SBG  ratio  is  10  and  4  respectively. 

Fig.  2.5.  Growth  of  large  grains  at  textured  band  boundaries,  as  predicted  by  MC 
modeling  for  BW/SBG  =  10  case,  after:  (a)  10  MCS;  (b)  20  MCS;  (c) 
50  MCS.  Slow  growth  of  subgrains  within  the  bands  also  occurs. 

Fig.  2.6.  Growth  of  grains  as  predicted  by  MC  modeling  for  BW/SBG  =  4  case, 
after:  (a)  30  MCS;  (b)  100  MCS;  (c)150  MCS. 

Fig.  2.7.  Two  bee  lattice  orientations  within  the  modeling  volume  and  major  GB 
orientation  within  the  modeling  volume  (most  probable  GBI)  (a),  and 
(b)  initial  microstructure  (cross-section)  for  modeling  volume. 

Fig.  2.8.  Modeling  domain  initial  texture  ODF;  (a)  (intensities  in  MU  per 
angular  unit)  and  (b)  ODF. 
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Fig.  2.9.  Effective  GB  mobility  as  function  of  both  scalar  GB  misorientation 
TTC  component  and  its  relation  with  modeling  volume  texture 
components  (A)  and  (B)  (schematically). 

Fig.  2.10.  MC  predicted  microstructures  after  (a)  75  MCS;  (b)  200  MCS;  c)  400 
MCS;  and  d)  800  MCS. 

Fig.  3.1.  Cubic  MD  modeling  domain  of  2300  atoms  with  two  lattices  in 
different  misorientations. 

Fig.  3.2.  Spherical  MD  modeling  domain  of  1650  atoms  with  two  lattices  in 
different  misorientations  and  inclinations. 

Fig.  3.3.  Simple  pair  interatomic  potential  and  STGB  geometry  assumed  for  test 
MD  run. 

Fig.  3.4.  Geometry  of  the  [100]  STGB,  projection  on  (001)  plane  (not  relaxed). 
Fig.  3.5.  [100]  STGB  model  (relaxed). 

Fig.  3.6.  Relative  STGB  [100]  energy  as  function  of  misorientation,  MD 
simulation 

Fig.  3.7.  Result  for  bee  STGB  [100]  obtained  in  [8]. 

Fig.  3.8.  Geometry  of  the  simulated  GB  (central  horizontal  plane)  with 

misorientation  (a)  (30,  90,  30)  and,  (b)  (30,  19,  30),  central  part  of  the 
modeling  volume. 

Fig.  3.9.  GB  energy  for  (30,  Theta,  30),  Theta  =  [0,  90]  boundary.  Positions  1 
and  2  correspond  to  atomic  structures  (a)  and  (b),  respectively,  Fig. 

3.8. 

Fig.  3.10.  Introducing  GBI  into  GB  energy  estimation  procedure. 

Fig.  3.11.  GB  energy  for  (30,  0,  30)  (0,0,0)  misorientation;  GBI  =  [0  to  70]  deg. 
(a)  all  data  and  (b)  -  points  for  0  and  1  deg  GBI  are  excluded. 

Fig.  3.12.  Relative  GB  energy  for  (45,  Theta,  45),  Theta  =  [0,  20]  boundary. 

Marked  positions  correspond  to  CLS  patterns  presented  on  Fig. 3. 14 

Fig.  3.13.  Relative  GB  energy  for  (45,  Theta,  45),  Theta  =  [0,  90]  boundary. 
Marked  position  corresponds  to  E9  CLS  GB. 

Fig. 3. 14.  GB  for  (45,  Theta,  45)  (0,  0,  0)  CLS  visualization:  Theta  -  (a)  34.6  (b) 
38.6  and  38,8  (c)  deg. 

Fig  3.15.  GB  energy  for  (45,  Theta,  45),  Theta  =  [38,  39]  boundary.  Position,  (c) 
corresponds  to  E9  CLS  GB. 

Fig.  3.16.  GB  energy  as  function  of  GBI.  GB  misorientation  (30,  30,  30)  (0,0,0) 

(a) .  Not  inclined  GB  position  is  shown  by  bold  arrow  on  the  GB 
energy  misorientation  dependence  (b). 

Fig.  3.17.  GB  energy  as  function  of  GBI.  GB  misorientation  (45,  0,  0)  (0,0,0), 
rotation  around  100  axis. 

Fig.3.18.  GB  inclination  for  (0,0,0)(20,20,20):  GBM  (a)  co  =  0  deg,  §  =20  deg; 

(b)  co  —  20  deg,  <|>  =  20  deg 
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Fig.  3.19.  Relative  GB  energy  for  GBI  (0,  0,  0)  (20,20,20).  Inverted  color 
scheme  used  for  points  on  lines  AB  and  CD. 

Fig.  3.20.  Relative  GB  energy  for  sections  A-B  (a)  and  C-D  (b)  (see  Fig.  2). 

Marked  position  (bold  arrows)  correspond  to  geometry,  for  which  GB 
normal  orientation  is  close  to  1 10  direction  of  (20,20,20)  oriented  part 
of  bicrystal. 

Fig. 3.21.  Relative  GB  energy  smoothed  (adjusted-averaging)  color  map  for  GBI 
(0,  0,  0)  (20,20,20),  <|>  e  [0,  90]  and  o  e  [0,  90]. 

Fig.  3.22.  Relative  GB  energy  3D  plot,  scale  divided  by  factor  2  is  used  for  best 
imaging 

Fig.  3.23.  GB  Energy  maps  (J/Unit  Cell)  as  function  of  GBI  for  series  of  GBM 
(0,0,0  -  x,y,z)  :  a)  x,y,z  =  1,3,7;  b)  6,15,55  c)  8,2,6  d)  25,20,60  e) 
40,15,10  f)  50,  0,  10 

Fig.  4.1  Initial  GB  geometry  for  3D  MC  and  MD  simulation  (a)  cross-section 
and  (b)  3D  view. 

Fig.  4.2.  GB  energy  as  function  of  GBI  (used  as  GB  effective  mobility  in  3D  MC 
runs).  GB  misorientation  (30,  30,  30)  (0,0,0). 

Fig.  4.3.  MC  simulated  GB  migration  in  bicrystal  (3D  half-loop  geometry)  with 
inclination  independent  GB  effective  mobility  (a)  and  locally  GB 
inclination  dependent  effective  mobility  (b). 

Fig.  4.4.  Cross  sections  are  parallel  to  the  001  planes  of  the  “external”  lattice. 

Fig.  4.5.  Effective  GB  mobility  (Fig.  1.33)  for  (0,  0,  0  20,  20,  20) 

misorientation  -  proportional  to  the  GB  relative  energy  GBI 
dependence 

Fig.  4.6.  Analysis  of  GB  faceting  in  3D  after  4000  MCS.  Reconstruction  by 
voxels  in  gray  light  (a)  and  colored  facets  (b). 

Fig.  5.1.  GB  energies  for  misorientations  of:  (a)  2  deg.;  (b)  3  deg.;  (c)  10  deg.; 

(d)  13  deg.;  (e)  29  deg  and  (f)  45  deg  presented  as  dependence  on 
Bunge  angles  variant  number  in  corresponding  dataset  (proportional  to 
first  rotation  angle,  for  GBI  (20,  20) 

Fig.  5.2.  Averaged  for  all  orientation  variants  GB  energy  as  function  of 

misorientation  (a)  GB  Inclination  =  (20,  20);  (b)  GB  Inclination  =  (30, 
30). 

Fig.  5.3.  Bicrystal  with  a  [1  1  0]  asymmetric  tilt  boundary  [15].  The 

misorientation  is  defined  as  the  angle  0  (=  38,9  deg)  between  the  [0  0 
1]  directions  of  two  single-crystal  grains,  and  the  inclination  is  defined 
as  the  angle  cp  between  the  boundary  plane  and  the  bisector  of  the 
misorientation 

Fig.  5.4.  Dependence  of  the  tilt  grain  boundary  energy  on  inclination  cp:  (a)  MD 
estimated,  lines  correspond  to  deep  spline  averaging  (dash  green)  and 
bin  (20)  averaging  (solid  blue)  and,  (b)  experimentally  measured  in 

[15]. 

Fig.  5.5.  Relative  GB  energy  for  GBI  (0,  0,  0)  (20,20,20).  (a)  Born-Mayer  map 
(b)  -  Lennard-Jones  map  and  (c)  Lennard-Jones  map  rescaled  for  4 
levels 


PAGE  44 

PAGE  44 

PAGE  45 

PAGE  45 

PAGE  46-47 

PAGE  48 

PAGE  48 

PAGE  49 

PAGE  50 
PAGE  50 

PAGE  52 

PAGE  55 

PAGE  56 

PAGE  56 

PAGE  57 

PAGE  58 


STCU 


PROJECT  P-262  (EOARD  068002)-  FINAL  PROJECT  REPORT 


PAGE  9 


Fig.  5.6.  Relative  GB  energy  for  GBI  (0,  0,  0)  (20,20,20)  calculated  with 

Lennard-Jones  (solid  line)  and  Born-Mayer  (dashed  line)  potentials. 

Fig.  6.1.  Optical  (left)  and  EBSD  (right)  images  of  a)  solid  solutioned;  b)  cold 
drown;  c)  50%  Recrystalized  and  d)  100%  recrystalized  material. 

Fig.  6.2.  Predicted  recrystallization  kinetics  for  Cases  A  and  B. 

Fig.  6.3.  Simulated  microstructures  for  Cases  1  (left)  and  2  (right).  GB  color  is 
proportional  (red  GBs  >15  deg  mis,  other  <15,  blue  <  2  deg) 

Fig.  6.4.  LCB  95%  Rex.  EBSD  Grain  Boundary  Map  (400*500  pm)  shows  grain 
arrangements  (marked  with  green  color)  which  looks  similar  to  Case  4 
simulated  patterns  (right). 

Fig.  6.5.  Initial  (as-deformed)  microstructure. 

Fig.  6.6.  Dependence  of  GB  energy  on  GBI  for  high  symmetry  GBMs:  (a) 
(45?45?45)  and  (b)  (45,45,0). 

Fig.  6.7.  Dependence  of  relative  GB  energy  on  GBI  for  general  GBMs  from 
table  6.1 

Fig.  6.8.  Recrystallization  kinetics  for  Cases  1  and  2. 

Fig.  6.9.  Microstructural  evolution:  (a)  Case  1;  and  (b)  Case  2. 

Fig.  6.10.  Some  favorite  GB  planes  for  Case  2,  60  MCS. 

Fig.  6.1 1.  MC  Potts  Rex.  and  GG  software:  (a)  modeling  volume  analysis;  (b) 
searching  for  GB  energy  value  within  included  database. 

Fig  7.1.  Steady  state  GB  migration  (half-loop)  geometry  as  proposed  for  MD 
estimation  of  GB  mobility  (a)  and  (b)  -  geometry,  often  used  for 
experimental  GB  mobility  estimation. 

Fig.  7.2.  Experimental  GB  mobility  estimation  for  A1  alloy  [1 1].  (a)  low 

temperature,  essential  GBI  dependence  (three  nearly  flat  segments),  (b) 
high  temperature,  almost  no  GBI  dependence 

Fig.  7.3.  Example  of  spatial  orientation  of  the  rolled  sample  with  respect  to  the 
magnetic  field  H,  provides  essential  driving  for  GBs  force  due  to  the 
anisotropy  of  the  magnetic  susceptibility  of  titanium  (magnetic 
annealing  of  titanium,  [12]) 

Fig.  7.4.  Relative  GB  energy  for  (0,0,0)  (20,20,20)  misoriented  GB  as  function 
of  GBI  component  (|>.  Open  circles  are  GB  energy  values  when 
additional  driving  force  is  effective. 

Fig.  7.5.  Example  of  flat  relaxed  GB  in  MD  simulation 

Fig.  7.6.  Migrating  hemispheric  GB  geometry  used  in  MD  simulation. 

Fig.  8.1.  CLS  visualization  tool  (a)  and  CLS  patterns  changing  from  LABs  to 

HABs  (twist  GB,  0=  2  deg.  (b),  4  deg.  (c),  6  deg.  (e)  12  deg  and  (f)15 
deg). 
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INTRODUCTION 

The  final  microstructures  and  textures  of  most  technical  alloys  are  the  result  of  their  multistage 
thermomechanical  processing.  Hence,  characterizing  and  understanding  of  micro  structure  and  texture 
evolution  during  various  thermomechanical  processes  are  important  for  improving  their  behavior  and 
properties. 

Recently,  computer  modeling  of  micro  structural  evolution  has  established  itself  as  an  essential  tool 
for  predicting  final  (post  processing)  microstructures  and  properties.  The  most  important  components  of 
polycrystalline  materials  are  internal  interfaces  and,  especially,  grain  boundaries  separating  adjacent 
crystals  (grains)  with  identical  chemical  composition  and  crystallography,  that  differ  only  by  mutual 
misorientation  of  their  crystal  lattices.  The  grain  boundaries  in  polycrystalline  create  a  three-dimensional 
net  spreading  throughout  the  material  (the  main  feature  of  the  microstructure);  they  have  different 
properties  and  therefore,  they  substantially  contribute  to  the  behavior  of  the  material  through  forming 
material  grain  structure  and  crystallographic  texture.  Due  to  a  large  variety  of  possible  misorientations  of 
the  adjacent  grains,  as  well  as  of  the  orientation  of  the  actual  joining  planes,  there  is  a  huge  number  of 
grain  boundaries  possessing  different  structures  and  thus,  different  properties. 

A  weak  point  of  grain  growth  and  recrystallization  modeling  is  lack  of  reliable  information  on 
characteristics  of  grain  boundaries  (GB  energy  and  mobility),  which  are  critical  for  quantitative  modeling. 
It  was  assumed  in  all  previous  computer  models  (and  in  experiments  on  grain  boundary  migration  in 
bicrystals)  that  the  spatial  orientation  of  a  boundary  does  not  seriously  affect  GB  energy  and  mobility.  This 
is  not  true,  because  in  the  case  of  a  curved  GB  its  structure  (and  properties)  changes  along  the  boundary. 

Hence,  the  aim  of  the  project  was  to  build  a  really  predictive  model  of  recrystallization  and  grain 
growth.  A  number  of  problems  were  solved,  including  accounting  of  an  exact  grain  boundary 
crystallogeometry,  i.e.  local  GBI  together  with  three-parameter  GBM  as  arguments  for  local  GB  energy 
and  mobility,  development  of  modeling  procedure.  MD  simulation  together  with  recent  experimental  data 
analysis  were  applied  for  obtaining  effective  values  for  GB  energy  and  mobility  for  all  possible 
combinations  of  exact  GBM  and  GBI.  An  attempt  of  estimating  realistic  data  for  MC  modeling  concerning 
deformed  technical  Ti-based  alloy  (LCB)  (grain  and  sub  grain  structure,  GBM  and  GBI  distributions)  and 
introducing  them  as  starting  conditions  for  MC  modeling  has  been  performed. 


STCU  PROJECT  P-262  (EOARD  068002)-  FINAL  PROJECT  REPORT  FPR  PAGE  1 1 

1.  EXACT  GB  GEOMETRY. 

Stage  1  of  project  was  aimed  at  implementing  of  an  exact  description  of  grain  boundary  geometry 
into  the  MC  Potts  modeling  procedure  instead  of  commonly  used  scalar  GBM.  It  was  shown  how  two- 
parameter  grain  boundary  inclination  (GBI)  could  be  accounted  in  modeling  procedure.  GBMs  where 
presented  in  their  exact  matrix  form. 

In  all  previous  computer  models  of  microstructure  evolution  it  was  assumed  that  the  spatial 
orientation  of  a  boundary  does  not  seriously  affect  grain  boundary  energy  and  mobility.  However,  structure 
and  properties  of  a  curved  grain  boundary  changes  along  the  boundary.  Hence,  to  build  a  really  predictive 
model  of  recrystallization  and  grain  growth  not  only  three-parameter  grain  boundary  misorientation 
(GBM)  but  local  grain  boundary  inclination  (GBI)  must  be  accounted  as  arguments  for  local  grain 
boundary  energy  and  mobility. 

In  this  project  formulation  of  the  mathematical  background  for  incorporating  the  exact  grain 
boundary  crystalogeometry  into  the  computer  simulation  routines  (algorithm  aimed  to  derive  local  grain 
boundary  inclination  (GBI)  for  the  Monte-Carlo  (MC)  Potts  modeling  volume  unit  (MU)  at  grain  boundary 
(GB))  was  carried  out.  Simple  and  effective  routine  for  GBI  estimation  was  created  and  implemented  into 
the  basic  MC  Potts  code.  GBM  was  implemented  into  MC  Potts  modeling  in  matrix  form.  Test  programs 
were  run  to  show  how  the  “pure”  GB  properties  dependence  on  GBI  affects  both  the  isolated  grain 
shrinkage  and  the  modeling  microstructure  evolution. 

1.1.  Incorporating  grain  boundary  inclination  (GBI) 

Algorithm  aimed  at  deriving  local  GBI  for  the  MC  Potts  modeling  volume  point  at  GB  was 
formulated  basing  on  following  assumptions.  Unlike  GBM,  the  GBI  should  be  determined  locally.  The 
GBI  is  a  subject  to  constant  change  while  GB  migrates,  hence  information  obtained  during  previous  MC 
steps  cannot  be  used  for  next  ones.  Therefore,  it  is  not  worthwhile  to  use  the  most  precise  approach, 
namely  GB  approximation  by  spatial  (3D)  analytical  function  because  one  needs  to  repeat  this  time- 
consuming  procedure  after  each  elementary  MC  event.  Another  problem  with  GBI  estimation  procedure, 
if  based  on  analytical  function  fit,  is  that  GB  may  consist  from  very  different  number  of  MUs;  up  to 
thousands;  this  makes  developing  of  the  unified  fit  procedure  very  complicated.  Let  us  consider  the  case 
when  the  GB  is  represented  by  function  z  =  f(x,y).  We  suggest  to  estimate  GBI  locally  using  planar 
linear  approximations  for  basic  Cartesian  GB  cross-sections.  Cross-sections  parallel  to  xz  andyz  through 
the  M0(x,  y,  z)  point  in  the  modeling  volume  for  which  local  GBI  should  be  estimated  are  shown  in  Fig. 
1.1. 

Partial  derivatives  df/dx  and  df/dy  define  tangential  lines  for  each  cross-section.  Then,  tangential 
plane  for  z  =  f(x,y)  includes  two  tangential  lines  within  the  cross-sections  parallel  to  xz  and  yz,  angles  a 
and  P  are  determined  by  linear  approximation.  Tangential  plane  normal,  that  characterizes  GBI  in  point 
M0  is: 


rn  = 


and  normal  line  (which  includes  given  MU)  equation  is: 
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i  —  Jq  V  —  Va  z  —  Zq, 

Next  algorithm  was  developed.  For  each  cross-section  the  nearest  vicinity  (48MUs)  of  the  MU  under 
consideration  was  taken  into  account,  as  shown  in  Fig.  1.2a.  Although  two  Cartesian  cross-sections  are 
sufficient  for  linear  approximation  sometimes  number  of  MUs  within  one  of  them  is  not  enough,  as  shown 
in  Fig  1.2b.  We  propose  to  use  three  Cartesian  cross-sections  to  guarantee  that  at  least  two  of  them  will  be 
suitable  for  linear  approximation. 


Fig.  1.1.  Geometry  of  two  Cartesian  cross  sections,  namely  xz  andyz,  for  the  GB  z  =J[x,y)  with  GBI  close  to 
the  xy  plane  . 


Aa 

f 

AVJ-Q 

-3-2-10  1  2  3 


b) 


Fig.  1.2.  GB  cross-sections:  (a)  cross  section  showing  MUs  neighboring  the  M0,  the  best  linear  fit  with  the 
normal  direction  projection;  and  (b)  cross  section  including  single  MU. 


For  procedure  error  estimation  a  model  of  spherical  grain  of  different  radius  R  was  tested.  Such 
test  geometry  insured  that  all  possible  GBIs  and  curvatures  would  be  tried.  The  second  advantage  is  that 
unlike  for  real  GB,  exact  GB  tangential  plane  position  is  known. 
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a)  correct  error 

c) 

Fig.  1.3.  Cross-sections  generated  using  test  model  which  represents  accuracy  of  the  GBI  estimation:  (a) 
exact  result  for  R=50MU;  (b),  (c)  error  origin  for  large  and  small  GB  curvature  radius,  respectively. 

It  was  assumed  that  acceptable  average  error,  i.e.  angle  between  exact  and  estimated  normals  is  2  deg., 
i.e.  not  more  than  commonly  used  Euler  angle  step  for  discrete  representation  of  orientations.  Fig.  1.3a 
presents  the  cross-section  of  simulated  spherical  grain  and  a  normal  line  projection  estimated  by  using  the 
above  procedure. 

The  acceptable  errors  were  obtained  for  more  than  99%  estimations  at  all  curvatures  having  radius 
above  20MU.  For  small  curvatures  a  number  of  successful  estimations  dropped  to  94%.  The  average  errors 
are  shown  in  Table  1.1.  For  GB  curvatures  above  50MU  average  error  in  GBI  estimation  remains  constant, 
i.e.  result  is  not  improved  with  further  growth  of  GB  curvature.  Figures  1.3b  and  1.3c  illustrates  origin  of 

errors  in  GBI  estimation  for  large  and  small  GB  curvatures,  correspondingly. 

1  . 


Fig.  1.4.  GB  curvature  radius  for  most  microstructures  is  larger  than  average  grain  size. 

Usually  GB  curvature  radius  20MU  corresponds  to  micro  structure  of  average  grain  size  10-12MU  [5]. 
Fig.  1.4  illustrates  correspondence  between  GB  curvature  in  real  modeling  microstructure  and  the  curvature 
of  the  spherical  grain.  The  conclusion  is  that  procedure  developed  cannot  define  correctly  GBI  for  small 
grains  (GB  segments)  because  the  whole  GB  segment  is  smaller  than  working  field  shown  in  Fig.  1.2a.  For 
fine  microstructures  we  propose  to  use  5x5  working  field  (25  MUs),  what  still  gives  acceptable  results. 
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It  should  be  mentioned  that  contrary  to  study  [7]  in  which  method  of  estimation  of  GBI  in  discrete 
micro  structural  modeling  for  modeling  volumes  was  proposed,  our  procedure  does  not  require  adding  extra 
information  for  data  storage. 


Table  1.1.  Error  (angle  between  estimated  and  exact  normal  to  GB)  statistics  versus  GB  curvature  radius  R, 
based  on  200  estimations  for  randomly  selected  GBI  for  each  R. 


Curvature  Radius,  R 

100 

50 

30 

25 

20 

15 

10 

8 

Average  error, 
deg 

1.641 

1.687 

2.255 

2.536 

2.423 

4.301 

6.562 

7.24 

Standard  deviation, 
deg 

1.029 

1.039 

1.642 

1.77 

1.643 

3.002 

4.301 

4.167 

Fraction  of  exact 

GBI  determination,  % 

10.5 

13.0 

11 

11 

10 

8 

5 

6 

Unsuccessful ,  % 
(error  >  2%) 

N/A 

N/A 

0.5 

0.5 

1.0 

0.5 

4.0 

6.0 
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1.2  Test  run  on  GBI  estimation  procedure 

For  procedure  test  run  model,  shrinkage  of  spherical  grain  (R  =  100MU)  was  modeled.  Such  test 
geometry  insured  that  all  possible  GBIs  and  curvatures  would  be  tried.  Test  mimics  the  one  described  in 
[4],  which  was  aimed  to  prove  the  MC  procedure  validity  by  modeling  the  shrinking  spherical  grain.  The 
grain/matrix  misorientation  was  taken  as  35  degrees.  In  [4]  the  GBM  was  the  only  parameter  which 
defined  the  boundary  mobility.  The  grain  remained  spherical  until  full  disappearing.  In  our  case,  the 
boundary  is  defined  not  only  by  GBM  but  GBI  as  well.  The  experimental  data  were  found  on  how  GBI 
affects  effective  GB  mobility/energy  [3].  Unlike  previous  works,  the  paper  [3]  was  not  restricted  to  plane 
boundaries  only  (simple  tilt  or  twist)  around  high  symmetry  lattice  directions.  Fig.  1.5  illustrates  these 
experimental  data  for  combined  twist  and  tilt  GBs  in  NiAl  alloy.  Our  test  MC  model  was  run  using  GB 
properties  versus  GBI  according  to  [3]. 


a: 

Pure  i 
Twist 


Pure  Tilt  0 

0 


30  45 

15  ^ 

Disorientation  angle  [deg] 


1  -  TTC 


c)  d) 

Fig.  1.5.  GB  effective  mobility  v.s.  GBI  [3]:  (a)  Optical  microscope  image  of  GB;  (b)  Dependence  of  the 
relative  GB  energy  on  crystallographic  parameters  of  the  GB  plane  in  two  adjacent  grains  (I  and  II). 
The  letters  A-H  refer  to  eight  positions  along  GB  and  (c)  Dependence  of  the  relative  GB  energy  on 
the  tilt/twist  component  (TTC)  and  disorientation  angle;  d)  averaged  effective  GB  mobility  as  TTC 
function  for  30-45  deg  misorientation  interval. 


For  convenience,  we  followed  the  tilt/twist  component  (TTC)  definition,  according  to  [3], 
TTC  =  \U,  V ,  W]  ’  n2 ,  where  both  vectors  were  normalized  to  unity.  Its  value  varied  from  0  for  pure  tilt  GBs 
to  1  for  pure  twist  GBs.  The  general  case  of  0  <  TTC  <  1  described  a  mixed  type  GB,  with  TTC 
representing  the  fraction  of  the  twist  component  in  the  GB.  GBM  between  grain  and  matrix  was  chosen  to 
be  between  30  and  45  deg  (this  value  to  be  referred  to  fig.  1.5c,  from  where  averaged  relative  GB  energy 
dependence  on  TTC  was  derived  (Fig.  1.5d)).  As  expected,  account  of  GBI  resulted  in  different  GB 
migration  velocities  for  different  TTC,  as  shown  in  Fig.  1.6. 
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a)  b)  c) 

Fig.  1.6.  Cross-sections  generated  using  test  model  which  represent  initially  spherical  grain  shrinkage  due 
to  relative  GB  energy  dependence  on  GBI  (TTC):  (a)  initial  state  R=100MU;  (b),  (c)  grain  shape 
after  50  and  100MCS.  Deviation  from  the  spherical  shape  is  evident. 

1.3  Implementation  of  the  GBI  estimation  procedures  into  3D  Monte 
Carlo  modeling  software 

Procedure  developed  for  GBI  estimation  cannot  define  correctly  GBI  for  smaller  grains  and  GB 
segments  close  to  triple  lines  because  the  whole  GB  segment  is  smaller  than  working  field.  Hence,  if 
M0  point  is  located  close  to  triple  line  or  triple  line  junction,  it  is  impossible  to  derive  GBI  from  the 
MC  Potts  modeling  volume  due  to  low  resolution  of  the  discrete  modeling  volume  grid.  Figure  1.7b 
illustrates  spatial  limitations  for  the  GBI  determination  procedure.  Since  the  procedure  for  GBI 
determination  needs  at  least  4  nearest  to  M0  GB  MUs  to  be  present  within  each  Cartesian  cross-section 
it  can  not  work  with  GB  curvature  radius  less  then  4  MU.  Areas  of  triple  lines  or  triple  line  junctions 
are  shown  in  Fig.  1.7b  by  deep  blue  filling.  The  procedure  was  proposed  that  excludes  triple  line 
vicinity  areas  from  the  GBI  estimation. 


a)  b) 

Fig.  1.7.  Modeling  volume  cross-sections:  (a)  cross  section  showing  GB  normal  estimated  at  M0  (projection 
within  the  section  plane;  and  (b)  the  same  cross-sections  with  several  more  GB  normal  lines 
projections,  at  which  areas  excluded  from  GB  normal  estimation  are  shown  by  deep  blue  color. 
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Grain  approximated 
with  the  Dodecahedron 


Grain  size,  MU 


b) 


Fig.  1.8.  (a)  Grain  shape  approximated  with  the  Dodecahedron  and  (b)  estimation  of  GB  volume  part  for 
which  GBI  can  be  derived. 


Excluded  part  of  GB  estimated  by  taking  the  Dodecahedron  as  the  grain  shape  approximation  and 
excluding  the  triple  line  volume  from  the  GB  volume  is  shown  in  Fig  1.8a.  Taking  into  account  that 
Dodecahedron  internal  radius  (1/2  of  grain  size)  is  approximately  1.12  S  (S  is  the  edge  of  polyhedron),  GB 
part  for  which  GBI  can  be  estimated  was  derived.  As  it  seen  from  Fig.  1.8b,  the  whole  algorithm  works 
good  starting  from  grain  sizes  of  20-25  MUs,  when  not  less  than  70%  of  GB  volume  for  most 
microstructures  can  be  accounted. 

Figure  1.9  presents  an  implementation  of  the  effective  GB  energy  which  depends  on  GBI  as  shown 
on  Fig.  1.5d.  Effect  of  GBI  on  microstructure  is  obvious,  however,  grain  shapes  look  unrealistic  due  to 
absence  of  misorientational  dependences  of  GB  properties.  Thus,  effect  of  “pure”  GBI  influence  on  GB 
kinetics  was  illustrated.  Modelling  with  accounting  of  both  GBM  and  GBI  dependencies  will  be  described 
in  section  1.5. 
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a) 


b) 


c) 


Fig.  1.9.  Cross-sections  generated  using  model  of  normal  grain  growth  with  only  GBI  dependence  of  GB 
effective  energy:  (a)  initial  microstructure;  (b,  c)  after  20  and  200  MCS  correspondingly. 


1.4.  Implementation  of  GBM  in  matrix  form 


The  MC  Potts  software  for  recrystallization  and  grain  growth  was  updated  with  the  GBM 
presentation  in  the  exact  matrix  form.  Previously,  the  single  parameter  consisting  of  the  lattice 
misorientation  between  neighboring  grains  was  used  to  define  the  relative  grain-boundary  energy  and 
mobility. 


s  =  arccos 


m , ,  +  777  2 2  +  m33  - 1 


V 


The  present  modeling  routine  uses  the  exact  matrix  misorientation  M  between  grains  having 
orientations  gi  and  g2 : 


M  =  M„  xM"1 

g\  g  2 


mn 

mn 

m 13' 

m2 1 

m22 

m23 

mn 

mn 

m33  J 

where  the  rotation  matrices  M  ,  Mg2  correspond  to  orientations  gj  and  g2  of  the  neighboring  grains.  If 

the  lattice  symmetry  is  higher  than  triclinic,  one  also  needs  to  take  into  account  the  existence  of  equivalent 
rotations. 

The  evaluation  of  the  matrix  misorientation  would  require  very  long  time  if  done  for  each  MC  trial, 
even  using  effective  calculation  schemes  such  as  that  based  on  quaternions  [5].  Therefore,  a  database  of  M 
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as  a  discrete  function  of  gj  and  g2  was  calculated  and  loaded  into  the  computer  memory  prior  to  modeling. 
Thus,  only  a  single  search  operation  should  be  performed  for  each  MC  trial.  Because  of  lack  of  reliable 
information  on  characteristics  of  grain  boundaries  (misorientation  dependent  GB  effective  mobility),  for 
present  test  simulations  the  schematic  matrix  form  of  GBM  dependence  has  been  developed,  which  is 
equivalent  to  scalar  form  of  GBM  dependence  presented  on  Fig.  1.10b.  Taking  into  account  that  scalar 
misorientation  between  neighboring  grains  is: 


£  =  arccos 


mu  +m12  +  m33  -1 


v 


grain  boundaries  were  classified  as  LABs  if  |  Sp(M)  |  =  I  mii+m22+m33 1  >  2.96,  transition  boundaries 
(from  LAB  to  HAB)  if  2.92<  |  Sp(M)  |  <  2.96  and  HABs  if  |  Sp(M)  |  <  2.92. 


As  it  is  shown  in  Chapter  2,  GBM  dependence  can  be  derived  from  molecular  dynamics  modeling. 


a)  b) 

Fig.  1.10.  Effective  GB  mobility  as  function  of  TTC  component  (a)  and  of  scalar  GB  misorientation  (b). 

Figure  1.11  (a,  b,  c)  presents  evolution  of  microstructure  modeled  with  the  above  assumption.  As 
expected,  more  realistic  final  microstructures  were  obtained  when  both  GBM  and  GBI  dependencies  were 
introduced  to  the  MC  routine  in  comparison  with  the  case  when  only  GBI  dependence  of  effective  GB 
mobility  was  accounted.  Some  lower  grain  growth  rate  was  observed  due  to  the  decrease  of  the  average 
GB  effective  mobility  when  GBM  is  accounted. 

If  initial  microstructure  (Fig.  1.11a)  evolves  only  under  GBI  dependence,  it  changes  as  shown  in  Fig. 
1.11  d,  e. 
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Fig.  1.11.  Cross-sections  generated  using  model  of  normal  grain  growth:  (a)  initial  micro  structure;  (b)  with 
no  misorientational  but  only  GBI  dependence  of  GB  effective  energy  after  200MCS;  and  (c,d)  with 
both  misorientational  and  GBI  dependence  of  GB  effective  energy  after  200  and  800  MCS, 
correspondingly. 


b)  b)  c) 

Fig.  1.12.  Effective  GB  mobility  as  function  of  scalar  GB  misorientation  only  (a);  of  TTC  component 
(b)  and  both  misorientation  and  TTC  component  (c). 

1.5.  MC  modeling  of  Grain  Growth  with  implemented  GBM/GBI 


To  clarify  the  GBI  effect  on  grain  growth  kinetics  initial  microstructure  for  three  comparative  cases  (A 
-  C)  was  taken  as  shown  on  Figure  1.13a.  :  Effective  GB  mobility  for  cases  A  to  C  was  defined  as: 

Case  A  -  GB  mobility  as  function  of  GB  misorientation  only,  no  GBI  influence  was  assumed  -  GBM 
usual  assumption  was  used  [4],  i.e.  low  angle  boundaries  were  assumed  to  perform  low  mobility  in 
comparison  with  the  high  angle  boundaries  (Fig.  1.12a); 

Case  B  -  GB  mobility  M=MGbm*MGbi  as  function  of  both  GBM  and  GBI  depending  on  tilt-twist 
component  (TTC)  averaged  for  35-45  deg  GB  misorientation  interval  (Figure  1.12b). 

Case  C  -  the  same  as  for  the  Case  B,  but  GBI  dependence  of  MGBi  was  approximated  with  smooth 
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surface,  presented  in  Figure  1.12c,  based  on  experimental  data  [3]. 

Figure  1.13b  presents  grain  size  distribution  of  the  modeling  domain  assumed  for  Cases  A  to  C  prior  to 
grain  growth  start.  The  micro  structure  has  the  1000  grains  within  the  2503  MU  domain;  hence  average 
initial  grain  size  is  31MU.  For  nearly  equiaxed  grain  shape  dispersion  of  GBI  distribution  was  relatively 
wide. 

Modeling  volumes  were  not  textured  initially;  hence,  the  effect  of  GBM  on  GB  motion  for  Cases  A  to 
C  was  the  same. 

As  expected,  predicted  grain  growth  kinetics  for  the  Case  A  is  so  called  “normal  grain  growth” 
kinetics,  for  which  average  grain  size  is  proportional  to  the  square  root  of  time  (Fig.  1.14).  For  Cases  B 
and  C,  at  initial  stage  grain  growth  is  slower  than  for  Case  A,  due  to  the  lower  GB  effective  mobility. 


a) 


b) 


Fig.  1.13.  Initial  microstructure  (a)  and  initial  grain  size  distribution  (b)  for  Cases  A  to  C. 


Fig.  1.14.  Predicted  grain  growth  kinetics  for  Cases  A  to  C. 
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for  boundaries  with  GBI  varied  mostly  around  pure  tilt  direction.  Effect  is  more  pronounced  for  the 
Case  C.  After  about  100  MCS  grain  growth  became  faster  for  both  Cases  B  and  C  due  to  selection  of  grain 
boundaries  with  high  effective  mobilities  but  then  slowed  again.  Hence,  the  periodic  alternating  of  faster 
and  slower  grain  growth  can  be  expected  for  longer  times,  as  in  strongly  textured  material  [4].  It  can  also 
be  seen  that  grain  growth  kinetics  in  the  cases  B  and  C  are  less  “smooth”,  as  compared  to  the  case  A. 

Initial  microstructure  for  the  Case  D  was  formed  with  the  same  average  grain  size  as  for  cases  A  to  C, 
but  having  1:3:3  aspect  ratio  (Figure  1.15a),  hence  it  performed  a  relatively  narrow  GBI  dispersion,  as 
compared  to  the  equiaxed  microstructure  (Figure  1.16).  The  microstructure  was  formed  with  nucleation 
and  grain  growth  technique  which  produced  initial  grain  size  distribution  as  shown  in  Figure  4b.  It  takes 
essential  time  for  the  first  several  grains  to  vanish  for  the  grain  size  distribution  like  on  Fig.  3b.  As  the 
result,  no  growth  of  average  grain  size  was  observed  for  the  case  D  (Figure  1.17).  Other  result  was  that  for 
the  material,  with  grain  shape  aspect  ratio  other  than  1:1:1.  Due  to  specific  dispersion  of  GBI  distribution, 
kinetics  of  grain  growth  in  non  equiaxed  micro  structure  has  proved  to  be  much  slower  (Figure  1.17).  It 
was  observed  that  the  microstructure  gradually  transformed  into  equiaxed  (Figure  1.18).  It  can  be  expected 
that  after  reaching  1:1:1  aspect  ratio,  the  kinetics  of  grain  growth  for  both  types  of  initial  microstructure 
will  follow  the  similar  dependence,  although  the  delay  at  the  beginning  stage  of  the  growth  will  result  in 
smaller  average  grain  size.  It  should  be  also  mentioned  that  some  initial  delay  in  grain  growth  in  non- 
equiaxed  microstructure  was  due  to  initial  grain  size  distribution  in  which  small  grains  were  absent  and 
therefore,  it  took  some  time  for  first  grains  to  disappear. 


Fig.  1.1 5.  Initial  microstructure  (a)  and  grain  size  distribution  (b)  for  Case  D. 
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Random  Matrics 


Cases  A-C 


Case  D 


^  Pure  tilt  GBf 


GBI 

Dispersion 


ABC 


Fig.  1.16.  Schematically:  difference  in  GBI  dispersion  for  Cases  A  to  C  and  D. 


Fig.  1.17.  Predicted  grain  growth  kinetics  for  Case  D  in  comparison  with  the  Case  C  (dashed  line). 
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d) 


Fig.  1.18.  Predicted  microstractures  for  Case  D  (a)  initial,  after  150  MCS  (b)  and  500  MCS  (c).  (d,  e). 
Initial  (d)  and  predicted  grain  size  distribution  (500  MCS). 


1.6.  MC  Modeling  of  Recrystallization  with  implemented  GBM/GBI 

Special  modeling  volume  was  constructed  to  study  how  implementation  of  GBI  dependence  of 
effective  GB  mobility  affect  recrystallization. 

The  initial  as-deformed  microstructure  was  formed  with  1:3:3  aspect  ratio  (Figure  1.19)  thus 
providing  the  narrow  GBI  dispersion.  The  model  domain  had  a  sharp,  single  component  initial  texture 
strong  enough  to  obtain  most  grains  similarly  oriented.  The  oriented  continuous  nucleation  of 
recrystallized  grains  was  assumed.  New  grains  were  assumed  to  be  misoriented  30  to  45  degrees  to  the 
deformed  matrix,  but  rotation  axes  were  randomly  selected.  Furthermore,  it  was  assumed  that  grain 
boundaries  were  preferred  sites  for  nucleation  of  recrystallized  grains.  The  nucleation  rate  for  each  MCS 
within  the  unrecrystallized  material  was  set  as  10-4  nuclei  per  MU3.  The  recrystallization  nuclei  were 
assumed  to  be  equiaxed  with  an  average  volume  of  15  MU3. 

Initial  stored  energy  distribution  is  shown  on  Figure  1.20.  Stored  energy  was  assumed  to  be  dense 
close  to  GBs  descending  to  grain  interior.  Recrystallization  front  velocity  was  assumed  to  be  proportional 
to  the  local  stored  energy  density,  as  well  as  to  the  effective  mobility.  In  the  reference  case  A  the  GB 
mobility  was  assumed  as  function  of  GB  misorientation  only.  Hence,  GBM  usual  assumption  was  used, 
i.e.  low  angle  boundaries  (LABs)  were  assumed  to  perform  low  mobility  in  comparison  with  the  high 
angle  boundaries  (HABs).  In  the  case  B,  GB  mobility  was  defined  by  GBM  and  GBI  dependencies 
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approximated  with  smooth  surface  based  on  experimental  data  taken  from  Ref.  3;  (Figure  1.21). 


a) 


b) 


Fig.  1.1 9.  (a)  3D  image  of  selected  deformed  grains  and  some  rex.  nuclei,  and  (b)  cross-section  (b) 
generated  for  simulation  initial  micro  structure. 


Fig.  1.20.  Stored  energy  distribution  shown  with  color  (blue  to  red  scale)  within  the  modeling  domain. 
Cross-sections  generated  transverse  plane  (a)  and  rolling  plane  (b)  and  3D  image  (c)  showing 
some  rex.  nuclei  and  not  recrystallized  volume  after  50  MCS  (Case  A). 


Misorientation,  deg 


Fig.  1.21.  Effective  GB  mobility  as  function  of  both  scalar  GB  misorientation  TTC  component. 
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Fig.  1.22.  MC  predicted  microstructures  after  75  MCS:  (a)  Case  A  (90%  rex.)  and  (b)  Case  B  (80%  rex.) 

Microstructural  changes  during  recrystallization  for  cases  A  and  B  after  equal  number  of  MCS  are 
illustrated  with  Figure  1.22.  For  assumed  recrystallization  model  (oriented  nucleation  at  grain  boundaries 
and  stored  energy  descending  from  grain  boundary  to  the  centre  of  the  deformed  grain)  average  GBI  for 
the  recrystallization  front  is  about  the  same  as  for  the  as  deformed  material  (Figure  1.22b).  The 
recrystallization  kinetics  depends  on  recrystallization  front  effective  mobility.  Hence  it  can  be  expected  to 
should  be  slower  when  GBI  is  accounted.  In  fact,  it  was  shown,  that  taking  into  account  GBI  essentially 
changes  recrystallization  kinetics,  and,  in  particular,  time  for  recrystallization  to  complete  (Figure  1.23). 
Figure  1.24  shows  the  microstructures  for  Cases  A  and  B  after  recrystallization  completes.  It  is  clear  seen 
that  the  average  grain  size  for  case  B  is  about  two  times  bigger  than  for  case  A  (18  and  32MU 
correspondingly). 


Fig.  1.23.  Predicted  recrystallization  kinetics  for  Cases  A  and  B. 
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Fig.  1.24.  Predicted  microstructures  for  Cases  A  and  B  after  recrystallization  completes:  a)  Case  A,  100 
MCS  b)  Case  B,  250  MCS. 
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2.  MODELING  OF  INHOMOGENIOU SLY  TEXTURED  MATERIAL 


Interface  for  introducing  local  textures  into  the  modeling  microstructure  was  developed.  Test 
programs  were  created  and  run  to  show  how  the  “pure”  GB  properties  dependence  on  GBI  affects  the 
special  model  microstructure. 

Basic  software  was  the  same  as  described  in  [1].  Applications  described  below  were  designed  using 
Borland  Delphi  tools  to  run  under  MS  Windows  2000/XP  OS. 

2.1  Introducing  texture  inhomogeneities  into  3D  modelling  volume 

Interface  for  introducing  the  spatial  texture  inhomogeneities  was  developed.  Figure  2.1  illustrates  the 
first  compiled  interface  version;  it  is  now  possible  to  divide  the  whole  modeling  volume  on  to  the  number 
of  local  spatial  regions,  and  different  textures  can  be  assigned  for  each  spatial  region.  Certainly,  the  local 
textures  have  to  match  desired  integral  texture  and  GBM  distribution  function  in  the  modeling  volume. 


Fig.2.1.  Interface  for  introducing  the  spatial  texture  inhomogeneities  into  the  modeling  volume. 


Fig.  2.2.  Local  textures  implemented  into  the  model  of  deformed  material  microstructure.  Spatial  region 
corresponds  to  the  001  fiber  texture,  region  B  is  close  to  Oil  texture  component  and  region  C 
random  background  (undefined). 
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To  define  the  grain  orientation  with  respect  to  the  specimen  (or  model  domain)  axes,  three  Euler  angles 
(cpi,  O,  and  cp2  in  Bunge  notation)  were  used.  Specific  values  of  the  crystallographic  orientation 
g  =  {(px ,  O,  (p2 }  (0<  cpi  <  360,  0<  O  <  90,  0<  cp2  <  90)  were  associated  with  each  MU  and  presented  as  MU 
color  according  to  color  scheme  as  shown  on  Fig. 2. 2.  The  grain-boundary  position  was  associated  with  the 
space  between  two  sites  having  unlike  orientations.  Each  grain  was  characterized  by  a  volume  equal  to  the 
number  of  MUs  that  comprised  it.  Hence,  we  can  define  the  ‘grain’s  center  of  gravity”  and,  on  the  base  of  its 
position,  make  a  decision  which  spatial  region  this  particular  grain  belongs  to. 

A  special  procedure  was  used  to  create  the  initial  textural  state  within  each  spatial  region  of  the  modeling 
volume  (Fig.  2.2)  and  thus  to  replicate  the  specified  ODF  for  each  region  as  closely  as  possible.  This 
procedure  was  based  on  the  following  probabilistic  function: 


1  1 


’LviaMsiG,)) 

MD 


x(g(G,)) 


{i.  s(G,>n, 

[0, 


in  which  g(G.)  is  the  orientation  of  the  z'-th  MU,  V (G, )  is  the  number  of  MUs  oriented  as  g(Gz),  Ng 
is  the  total  number  of  orientations  g  within  the  volume  Qj,  and  VMD  is  volume  of  the  overall  model  domain. 
From  tests  on  model  textures,  it  was  determined  that  this  procedure  did  not  work  well  when  the  number  of 
grains  in  the  local  zone  was  less  then  50,  but  gave  very  good  results  when  there  were  100  or  more  grains. 


Interface  for  introducing  the  spatial  texture  inhomogeneities  was  used  for  creation  of  initial  modelir 
volumes  with  several  differently  textured  bands.  To  show  that  interface  together  with  corresponds 
program  code  can  be  used  to  simulate  the  role  of  scale  factors,  such  as  local  texture  band  width  (BW)  ( 
subgrain  size  (SBG),  in  microstructure  evolution,  volumes  with  two  essentially  different  subgrain  sizes  wei 
created. 

2.2.  MC  modeling  with  texture  inhomogeneities  in  3D  (accounting  GBM  only). 

Figure  2.3  illustrates  the  technique  of  initial  modeling  volume  creation  using  the  current  interfac 
version;  it  is  possible  to  divide  the  whole  modeling  volume  on  to  the  number  of  local  spatial  regions,  ar 
different  textures  can  be  assigned  for  each  spatial  region.  Adjusting  the  volumes  and  shapes  of  the  local] 
textured  areas  is  the  way  to  match  desired  integral  texture  and  GBM  distribution  function  in  the  modelir 
volume. 

Cross  sections  of  initial  modeling  volumes,  having  BW/SBG  ratio  10  and  4,  are  presented  on  figure  2.‘ 
Figures  2.5  and  2.6  presents  the  temporal  evolution  predicted  by  MC  simulation  for  both  cases,  the  GBI 
and  GBI  dependencies  were  the  same  as  used  for  normal  grain  growth  simulation. 

On  figures  2.3  to  2.6  different  colors/grayscale  levels  were  used  for  different  grain  orientations,  henc 
spatial  texturing  is  clearly  seen. 
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This  component  texture  will  be  placed  in  area  specified  by  structure  elements. 
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Fig.  2.3.  Local  textures  implemented  into  the  model  of  material  microstructure,  with  bands  consisting  from 
small  subgrains  of  low  misorientation. 


a) 


b) 


Fig.  2.4.  Cross  sections  of  locally  textured  initial  modeling  volumes:  (a)  and  (b)  BW/SBG  ratio  is  10  and  4 
respectively. 


Fig.  2.5.  Growth  of  large  grains  at  textured  band  boundaries,  as  predicted  by  MC  modeling  for 
BW/SBG  =  10  case,  after:  (a)  10  MCS;  (b)  20  MCS;  (c)  50  MCS.  Slow  growth  of  subgrains  within 
the  bands  also  occurs. 
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Fig.  2.6.  Growth  of  grains  as  predicted  by  MC  modeling  for  BW/SBG  =  4  case,  after:  (a)  30  MCS;  (b) 
100  MCS;  (c)150  MCS. 

Essential  differences  in  grain  growth  rate  and  final  microstructures  were  predicted  for  these  two  cases. 
Hence,  it  can  be  concluded  that  program  code  for  MC  simulation  of  recrystallization  and  grain  growth  is  a 
suitable  tool  for  studying  how  scale  factors,  such  as  local  texture  band  width  to  subgrain  size  ratio,  affect 
final  microstructures. 


2.3.  MC  Modeling  with  implemented  GBM/GBI  (Grain  Growth,  spatially 
inhomogeneous  texture) 

The  effect  of  grain  boundary  mobility  anisotropy,  caused  by  GBI,  on  grain  growth  kinetics  was 
shown  by  comparison  of  grain  structure  evolution  within  two  parts  of  the  modeling  volume  having 
different  sharp  one  component  textures  in  the  initial  condition.  For  such  modeling  volume  each  part 
initially  was  characterized  with  narrow  GBI  distributions  around  two  different  GB  inclinations. 

The  initial  microstructure  was  formed  with  1:5:5  aspect  ratio  (Figure  2.7b)  thus  providing  the 
narrow  GBI  dispersion,  around  the  plane,  indicated  on  Fig.  2.7a.  The  model  domain  had  a  sharp,  double 
component  initial  texture  (Fig.  2.8)  strong  enough  to  obtain  most  grains  similarly  oriented.  Microstructure 
of  modeling  domain  was  uniform,  but  texture  components  were  spatially  separated  and  form  locally 
textured  zones  (A)  and  (B).  GB  mobility  was  defined  by  GBM  and  GBI  dependence  approximated  with 
smooth  surface  based  on  experimental  data  taken  from  Ref.  3  (Figure  2.9).  Figure  2.9  also  demonstrates 
how  different  GBI  in  A  and  B  zones  result  in  different  effective  mobilities  for  nearly  fiat  LABs  within  the 
zones.  The  interface  between  A  and  B  zones  formed  by  HABs  with  relatively  high  effective  mobility. 

Predicted  evolution  of  the  modeling  volume  is  presented  on  Fig.  2.10.  It  is  obvious,  that  due  to 
lower  effective  mobility  within  the  B  zone,  grains  within  this  area  are  growing  slower  than  within  the  A 
zone.  Grains  aspect  ratio  decreases  with  time  in  both  zones  (but  slower  in  zone  B).  As  the  result,  the  GBI 
dispersion  increases,  i.e.  the  difference  between  effective  GB  mobilities  within  the  A  and  B  zones 
decreases.  After  200  MCS  both  zones  perform  grain  growth  kinetics  close  to  normal  grain  growth. 
However,  grains  within  the  A  zone  are  bigger;  as  a  result,  grains  of  zone  B  are  gradually  consumed  by 
zone  A  grains  and,  finally,  zone  B  vanishes. 

Hence,  it  was  demonstrated  essentially  different  grain  growth  kinetics  for  the  differently  textured 
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a) 

Fig. 2. 7.  Two  bee  lattice  orientations  within  the  modeling  volume  and  major  GB  orientation  within  the 
modeling  volume  (most  probable  GBI)  (a),  and  (b)  initial  microstructure  (cross-section)  for 
modeling  volume. 


Fig. 2. 8.  Modeling  domain  initial  texture  ODF;  (a)  (intensities  in  MU  per  angular  unit)  and  (b)  ODF 
projection  indicating  texture  components. 


Fig.2.9.  Effective  GB  mobility  as  function  of  both  scalar  GB  misorientation  TTC  component  and  its 
relation  with  modeling  volume  texture  components  (A)  and  (B)  (schematically). 
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c)  d) 

Fig.  2.10.  MC  predicted  microstructures  after  (a)  75  MCS;  (b)  200  MCS;  c)  400  MCS;  and  d)  800  MCS 
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3.  MD  ESTIMATION  OF  RELATIVE  GB  ENERGIES 

It  was  proposed  in  the  project  to  use  advanced  MD  simulation  as  a  source  of  the  data  on  grain 
boundary  energy.  A  MD  simulation  includes  selection  of  numerical  integration  scheme  and  employment 
of  appropriate  boundary  conditions.  The  modeling  domain  in  3D  with  periodic  boundary  conditions  and 
max.  total  number  of  atoms  of  about  105-106  was  utilized.  Usage  of  a  large  modeling  domain  makes  the 
cell  shape  restriction  less  important  compared  to  a  small  systems  with  a  few  hundred  atoms. 

Starting  from  the  simple  cases  of  tilt  and  rotation  boundaries  for  which  similar  calculations  have 
already  been  done  by  other  authors,  a  number  of  grain  boundaries  with  arbitrary  GBI  and  GBM  were 
studied. 

MD  applications  and  code  developed  were  designed  using  Borland  Delphi  tools  to  run  under  MS 
Windows  2000/XP  OS.  MD  simulation  utilizes  Born-Myer  and  Lennard-Jones  type  potentials  with  rigid 
boundary  conditions,  to  reproduce  differently  misoriented  bicrystal  with  different  GBIs.  GB  relaxation 
was  assumed  to  occur  in  two  stages  -  sliding  (lattice  translation  in  GB  plane  and  along  GB  plane  normal 
direction)  and,  finally,  relaxation  process  for  each  individual  atom. 

For  lattice  visualization  the  VMD  code  [6]  was  used  (as  a  freeware). 


3.1.  Testing  of  cubic  and  spherical  modeling  domains. 

First  runs  of  molecular  dynamic  routines  were  aimed  at  forming  of  bicrystals  of  desired  GMB  and  GBI 
in  both  cubic  and  spherical  modeling  domains.  As  the  first  step,  stable  lattices  were  obtained  with  using 
the  simple  Born-May er  potential.  Figure  3.1  illustrates  the  modeled  bicrystals  with  GB  misorientation  0, 
10  and  20  degrees  obtained  by  simple  tilt. 

Figure  3.2  shows  bicrystals  with  two  different  GBI  (0  and  90  deg  tilt)  and  modeled  by  using  the 
spherical  modeling  domain.  Utilizing  the  spherical  modeling  domain  insures  that  domain  boundary 
conditions  affect  GB  energy  n  the  same  manner  (not  depending  on  GBI)  -  what  is  essential  for  the  GB 
energy  calculating  routines. 


Fig.  3.1.  Cubic  MD  modeling  domain  of  2300  atoms  with  two  lattices  in  different  misorientations. 
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Fig.  3.2.  Spherical  MD  modeling  domain  of  1650  atoms  with  two  lattices  in  different  misorientations  and 
inclinations. 


First  test  run  of  the  MD  model  for  GB  energy  estimation  was  performed  using  3D  atomistic  model 
for  GB  (Fig.  3.2).  The  potential  energy  for  the  z-th  atom  was  calculated  using  Lennard-Jones  potential: 


4>i  =  4f‘  E(W^>12  -  W1'#)6}* 

where  8  and  a  are  the  material  parameters,  and  r^  is  the  distance  between  the  i-th  and  the  j-th  atoms.  The 
GB  bulk  energy  is  defined  as 


4>gb  = 

isV 


where  Sgb  is  the  area  of  the  GB,  and  the  summation  is  done  for  all  atoms  in  the  domain  V.  Domain  V  is 
chosen  around  the  GB,  excluding  the  outer  surfaces  of  the  modeling  volume. 

Next  parameters  were  used  for  test  simulation:  a  =  0.22  nm,  8  =  0.08  eV,  which  correspond  to  lattice 
parameter  a  =  0.28665  nm  and  elastic  modules  Cn  =  230  GPa,  Cu  =  135  GPa  and  C44  =117  GPa  [7]  to 
mimic  interatomic  potential  used  in  [8]  for  bee  iron  (Figure  3.3). 


Fig.  3.3.  Simple  pair  interatomic  potential  and  STGB  geometry  assumed  for  test  MD  run. 
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GB  energy  estimation  was  done  for  symmetric  tilt  [100]  GB  (STGB)  (Figure  3.4)  with  misorientation 
from  0  to  90  deg,  with  angular  step  of  2  deg.  GB  estimation  was  done  for  “relaxed”  grain  boundaries,  i.e. 
GBs  with  minimal  free  energy  for  given  crystallographic  parameters.  Relaxation  procedure  included  3 
stages:  a)  generating  two  lattices  in  given  misorientation  within  the  same  modeling  volume;  b)  rigid 
relaxation,  i.e.  mutual  shift  of  lattices  to  find  the  position  with  minimal  GB  energy,  and  c)  individual  atom 
relaxation  which  form  stable  GB  structure  (Figure  3.5).  Last  stage  was  realized  with  usual  MC  procedure 
-  randomly  selected  atom  was  allowed  to  move  for  small  finite  distance  towards  energetically  preferred 
position;  then  procedure  repeated  until  GB  energy  reached  its  minimum. 

Calculated  relaxed  STGB  energies  are  presented  on  Fig  3.6.  using  the  relative  values  between  0 
(monocrystal)  and  1 . 

As  it  seen,  calculated  energies  of  special  boundaries  are  generally  lower  what  can  accepted  as  a  proof 
the  adequacy  of  the  calculation  procedure,  developed  in  this  project.  The  results  for  this  particular  case 
well  correspond  to  the  results  of  Ref.  [9]  (Figures  3.6  and  3.7).  But  contrary  Ref.  [8]  and  other  similar 
publications,  our  approach  allows  to  calculate  energy  of  not  only  tilt  or  twist  boundaries  but  as  we 
suppose,  it  can  be  used  for  any  type  of  the  boundary. 


Fig.  3.4.  Geometry  of  the  [100]  STGB,  projection  on  (001)  plane  (not  relaxed). 


Fig.  3.5.  [100]  STGB  model  (relaxed). 


Additionally,  GB  energies  of  special  grain  boundaries  (Table  3.1)  were  calculated. 
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Table  3.1. 


GB  plane 

Misorientation,  deg 

X 

Minimum  relative  GB  energy 

210 

53.13 

5 

0.242 

310 

36.87 

5 

0.250 

320 

67.38 

13 

0.382 

510 

22.69 

13 

0.653 

530 

61.93 

17 

0.425 

410 

28.07 

17 

0.651 

Misorientation,  deg 


Fig.  3.6.  Relative  STGB  [100]  energy  as  function  of  misorientation,  MD  simulation. 


Fig.  3.7.  Result  for  bee  STGB  [100]  obtained  in  [8]. 


Good  correspondence  obtained  with  previous  results  for  bee  STGB  [100],  obtained  with  analytical 
embedded  atom  method  [8],  in  particular,  low  relative  energies  for  several  special  GBs. 
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3.2.  Grain  boundary  energy  as  function  of  GBM. 

Using  preliminary  results  on  MD  modeling  (see  above  3.1  section)  two  runs  of  the  MD  model  for  GB 
energy  estimation  were  performed  using  3D  spherical  atomistic  model  for  GB.  Significantly  bigger  than 
previously  modeling  volume  was  used  (2003  atoms).  The  first  run  was  the  case  of  generally  misoriented 
GB  and  the  second  one  was  to  estimate  GBI  dependence  for  GB  of  certain  misorientstion. 

GB  energy  estimation  was  done  for  general  non  symmetric  GB  (Figure  3.8)  with  misorientation  (30, 
0,  30)  when  0  is  changing  from  0  to  90  deg,  with  angular  step  of  1  deg  (another  grain  has  (0,  0,  0) 
orientation).  GB  estimation  was  done  for  “relaxed”  grain  boundaries,  i.e.  GBs  with  minimal  free  energy 
for  given  crystallographic  parameters. 

As  it  is  seen,  calculated  energies  of  special  boundaries  are  generally  lower  what  can  be  accepted  as  a 
proof  the  adequacy  of  the  calculation  procedure.  Although,  relatively  high  modeling  volume  was  used  in 
simulation,  there  was  no  correlation  between  the  coincide  lattice  period  E  and  minimums  depth  in  GB 
energy  (Figure  3.9).  This  might  be  because  of  the  rather  big  angular  step  of  1  deg,  which  special 
boundaries  are  not  correctly  defined. 


a)  b) 


Fig.  3.8.  Geometry  of  the  simulated  GB  (central  horizontal  plane)  with  misorientation  (a)  (30,  90,  30)  and, 
(b)  (30,  19,  30),  central  part  of  the  modeling  volume. 
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Fig.  3.9.  GB  energy  for  (30,  Theta,  30),  Theta  =  [0,  90]  boundary.  Positions  1  and  2  correspond  to  atomic 
structures  (a)  and  (b),  respectively,  Fig.  3.8. 
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3.3.  Grain  boundary  energy  as  function  of  GBI. 


m=> 

Relaxation 


Fig.  3.10.  Introducing  GBI  into  GB  energy  estimation  procedure. 

Proposed  3D  MD  approach  allows  calculating  the  GBI  dependence  of  GB  energy  at  certain  GBM. 
Schematically  the  routine  is  presented  on  Fig.  3.10.  As  an  example,  the  calculation  was  performed  for  bee 
lattice  with  (30,  0,  30)  (0,0,0)  misorientation. 

Much  higher  GB  energies  for  GBI  oriented  along  the  high  symmetry  lattice  directions  and  as 
compared  to  other  GBIs  was  demonstrated  (Figure  3.11a).  If  this  special  GBI  value  is  excluded,  the 
variation  of  GB  energy  for  a  given  GBM  was  found  to  be  not  wery  significant  (Figure  3.1  lb).  This  result 
is  in  good  correlation  with  data  on  effective  interplanar  spacing  vs  GBI  dependencies  [7,  9]. 
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Fig.  3.11.  GB  energy  for  (30,  0,  30)  (0,0,0)  misorientation;  GBI  =  [0  to  70]  deg.  (a)  all  data  and  (b) 
points  for  0  and  1  deg  GBI  are  excluded. 
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Fig.  3.12.  Relative  GB  energy  for  (45,  Theta,  45),  Theta  =  [0,  20]  boundary.  Marked  positions 
correspond  to  CLS  patterns  presented  on  Fig.3.14 

Figure  3.12  presents  GB  energy  estimation  for  general  type  type  non  symmetric  GB,  for  (45,  0,  45) 
(0,  0,  0)  misorientation  and  GBI  parallel  to  100  plane  of  bottom  part  of  bicrystal,  when  0  is  changing 
from  0  to  20  deg,  with  angular  step  of  1  deg.  It  is  clearly  seen,  by  comparison  with  corresponding  CLS 
patterns  (Fig.  3.14)  that  transition  from  LABs  to  HABs  occurs  through  decreasing  an  area  of  nearly 
coherent  patterns.  This  results  in  more  dislocations  and,  finally,  in  forming  HAB  within  the  relaxed  GB. 

Full  symmetric  twist  sequence  of  this  GBs  family  is  presented  on  Figure  3.13. 


Fig.  3.13.  Relative  GB  energy  for  (45,  Theta,  45),  Theta  =  [0,  90]  boundary.  Marked  position  corresponds 
to  E9  CLS  GB. 
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(a)  (b)  (c) 

Fig.3.14.  GB  for  (45,  Theta,  45)  (0,  0,  0)  CLS  visualization:  Theta  =  (a)  34.6  (b)  38.6  and  38,8  (c) 

deg. 


Misorientation,  deg 


Fig  3.15.  GB  energy  for  (45,  Theta,  45),  Theta  =  [38,  39]  boundary.  Position,  (c)  corresponds  to  £9 
CLS  GB. 

As  it  is  seen  from  Figures  3.13  and  3.15,  according  to  MD  simulations,  the  closest  to  £9  CLS  grain 
boundary  lays  between  38.7  and  38.8  degrees.  It  should  be  noted  here,  than  when  we  were  trying  to 
examine  smaller  angular  steps,  no  father  details  were  seen  on  graph  (Fig.  3.15),  obviously  due  to  the  finite 
size  of  the  modeling  volume. 

Additionally,  it  was  demonstrated  how  relative  GB  energy  is  changing  with  variation  of  GBI 
(Figure  3.16).  It  was  concluded  from  results  presented  on  Figure  3.16  that  almost  all  inclined  GBs  had  3 
to  7  times  less  energies,  that  not  inclined  ones,  i.e.  inclined  GBs  are  “better  relaxed”.  However,  difference 
in  GB  energy  for  different  inclinations  isn’t  so  essential.  It  was  shown  also  that  in  the  case  of  special 
misorientations  like  that  presented  on  Fig.  3.17,  inclination  does  not  change  the  boundary  energy  as 
compared  to  that  noticed. 
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Fig.  3.16.  GB  energy  as  function  of  GBI.  GB  misorientation  (30,  30,  30)  (0,0,0)  (a).  Not  inclined  GB 
position  is  shown  by  bold  arrow  on  the  GB  energy  misorientation  dependence  (b). 
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Fig.  3.17.  GB  energy  as  function  of  GBI.  GB  misorientation  (45,  0,  0)  (0,0,0),  rotation  around  100  axis. 

3.4.  Relative  GB  energy  dependencies  on  GBI  and  GBM.  2D  maps. 


Using  advanced  MD  modeling  software  (see  above  sections  3. 1-3.3)  investigation  of  GB  energies  as 
functions  of  GBM  and  GBI  was  being  continued  for  bee  bicrystal  for  general  case  of  not  special  GB 
geometry.  For  further  MD  and  MC  simulations  unified  notification  considered  introduced,  based  on  the 
Cartesian  axis,  coincided  with  100,  010,  and  001  directions  in  one  the  neighboring  grains  (the  only 
exception  will  be  global  texture,  based  on  the  sample’s  or  modeling  domain  coordinate  system).  GBI  was 
defined  with  GB  normal’s  angles  oo  and  4>  (see  Fig.3.18). 

As  starting  point,  GB  energy  was  calculated  for  (0,  0,  0)  (20,  20,  20)  GBM  and  GBI  parallel  to  100 
plane  of  part  of  bicrystal  with  orientation  (0,0,0). 

Using  above  misorientation,  map  for  relative  grain  boundary  energy  as  function  of  GBI,  (normal  to 
GBI  direction  defined  with  §  and  co  angles)  using  angular  step  of  2  deg  both  for  §  and  oo  (about  2000 
points)  was  built  (Fig.3.19). 

Figure  3.19  presents  relative  GB  energy  estimation  for  general  type  HAB,  for  4>  e  [0,  90]  and  oo  e  [0, 
90]  angular  intervals.  For  comparison  of  2D  map  with  actual  data  relative  GB  energy  along  A-B  and  C-D 


(a)  (b) 

Fig.3.18  GB  inclination  for  (0,0,0)(20,20,20):  GBM  (a)  co  =  0  deg,  4>  =20  deg;  (b)  oo  =  20  deg,  4>  =  20  deg. 
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Fig.  3.19.  Relative  GB  energy  for  GBI  (0,  0,  0)  (20,20,20).  Inverted  color  scheme  used  for  points  on 
lines  AB  and  CD. 


(a)  (b) 

Fig.  3.20.  Relative  GB  energy  for  sections  A-B  (a)  and  C-D  (b)  (see  Fig.  2).  Marked  position  (bold 
arrows)  correspond  to  geometry,  for  which  GB  normal  orientation  is  close  to  110  direction  of 
(20,20,20)  oriented  part  of  bicrystal. 


Figure  3.21  presents  actual  data  points  (Fig.  3.19)  as  slightly  smoothed  (adjusted-averaging)  color 
map.  As  it  is  well  seen  from  the  2D  map  -  peaks  for  GB  energy  corresponds  (are  close  to)  positions  of 
low-indexed  directions.  Hence,  it  can  be  concluded,  that  GB  energy  has  its  maximum  values  when  GBI 
normals  are  within  the  certain  crystallographic  zone  of  one  of  the  bicrystals. 
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To  form  full  database  of  GB  as  function  of  GBM+GBI  with  misorientation  step  of  2  degrees  for  one 
particular  metal  or  alloy  we  need  about  45*45*28  «  58  000  of  such  maps  to  be  generated  (See  Fig.  3.23). 
However,  when  result  is  verified  for  few  of  them,  calculations  can  be  completed  automatically.  Then, 
when  both  databases  for  relative  GB  energy  and  mobility  will  be  created,  it  can  be  used  in  3D  MC 
simulation  of  recrystallization  or  grain  growth  without  any  previously  used  simplifications,  such  as  Tilt- 
Twist  Component  (Fig  1.21). 
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Fig.3.21.  Relative  GB  energy  smoothed  (adjusted-averaging)  color  map  for  GBI  (0,  0,  0)  (20,20,20), 
(|>  e  [0,  90]  and  oo  e  [0,  90], 


Fig.  3.22.  Relative  GB  energy  3D  plot,  scale  divided  by  factor  2  is  used  for  best  imaging. 
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e) 


f) 


Fig.  3.23  GB  Energy  maps  (J/Unit  Cell)  as  function  of  GBI  for  series  of  GBM  (0,0,0  -  x,y,z) 
a)  x,y,z  =  1,3,7;  b)  6,15,55  c)  8,2,6  d)  25,20,60  e)  40,15,10  f)  50,  0,  10 
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4.  MC  MODELING  WITH  IMPLEMENTED  GBI  AND  MD  DATA  (STEADY 
STATE  GB  MIGRATION) 

Steady  state  GB  migration  (GB  misorientation  (30,  30,  30)  (0,0,0))  for  half-loop  GB  geometry 
(Figure  4.1)  was  studied  in  3D  with  MC  approach  for  two  cases: 

a)  inclination  independent  GB  effective  mobility  and 

b)  locally  GB  inclination  dependent  effective  mobility. 

Effective  mobility  was  assumed  to  be  proportional  to  relative  GB  energy  as  determined  early  with  MD 
approach  (Figure  4.2). 


(a) 


(b) 


Fig.  4.1  Initial  GB  geometry  for  3D  MC  and  MD  simulation  (a)  cross-section  and  (b)  3D  view. 
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Fig.  4.2.  GB  energy  as  function  of  GBI  (used  as  GB  effective  mobility  in  3D  MC  runs).  GB  misorientation 
(30,  30,  30)  (0,0,0). 

It  was  shown  that  GBI  dependent  effective  GB  mobility  assumed  on  the  base  of  MD  GB  energy 
calculations  when  applied  to  MC  simulated  GB  migration  results  in  GB  faceting  in  3D  (Figure  4.3).  Some 
tendency  for  GB  with  GBI  dependent  effective  mobility  faceting  can  be  already  seen  after  the  10th  MCS. 
Faceting  becomes  more  pronounced  during  GB  migration,  although  the  number  of  facets  decreases  with 
time.  No  evidence  of  faceting  was  found  for  the  case  (a). 
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Fig.  4.3.  MC  simulated  GB  migration  in  bicrystal  (3D  half-loop  geometry)  with  inclination  independent 
GB  effective  mobility  (a)  and  locally  GB  inclination  dependent  effective  mobility  (b). 
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Three  dimensional  analysis  of  GB  shape  evolution  from  the  initial  one  to  the  steady  state  migration 
condition  was  performed.  Initial  GB  configuration  0  MCS  (light-grey),  Figure.  4.4.  Map  for  effective  GB 
mobility  (Fig.  4.5)  for  (0,  0,  0  20,  20,  20)  misorientation  is  proportional  to  the  GB  relative  energy  GBI 
dependence.  GBI  normal/plane  (0,  0)  is  parallel  to  001  direction  for  “external”  lattice  (vertical 
direction/horizontal  plane)  on  Fig.  4.4a). 

GB  configuration  after  4000  MCS  (dense-grey)  and  orthogonal  cross  sections  analyzed, 
schematically  at  Table  4.1  and  Figure  4.6. 


(a)  (b) 

Fig.  4.4.  Cross  sections  are  parallel  to  the  001  planes  of  the  “external”  lattice. 


Inclination  angle  (f> 

Fig.  4.5.  Effective  GB  mobility  (Fig.  1.33)  for  (0,  0,  0  20,  20,  20)  misorientation  -  proportional  to  the  GB 
relative  energy  GBI  dependence. 
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Table  4.1.  The  sequences  of  migrating  GB  orthogonal  cross  sections. 
White  lines  represent  initial  configuration. 


MCS 
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Fig.  4.6.  Analysis  of  GB  faceting  in  3D  after  4000  MCS.  Reconstruction  by  voxels  in  gray  light  (a) 
and  colored  facets  (b). 
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5.  FORMATION  AND  ANALYSIS  OF  GB  PROPERTIES  DATABASE 


Database  for  GB  energy  as  function  of  exact  GB  geometry  was  created,  partially  filled  with  MD  data 
and  integrated  to  the  MC  Potts  modeling  software. 

Database  was  developed  in  a  form  of  electronic  table  for  (3-Ti  GB  properties  data,  namely  GB 
energy  and  GB  relative  energy.  Energy  values  are  presented  as  function  of  three  GB  misorientation  and 
two  GB  inclination  parameters  and  allow  data  filling  with  angular  step  1  degree. 

Structure  of  GB  energy  database  and  calculated  values  ranges  up  to  date  are  presented  on 
Table  5.1.  Angular  step  is  given  in  brackets  next  to  the  angular  range. 

“Calculated”  value  -  indicates  that  position  in  database  filled  by  MD  simulation  procedure, 
developed  in  the  frames  of  the  project.  However,  available  computational  resources  were  not  enough  to 
complete  all  database  -  not  yet  completed  GBs  marked  as  “In  Progress” 


Table  2.1. 


GBM  GBM 

Euler  (1)  Euler  (2) 


GMB 
Euler  (3) 


GBI(l)  GBI  (2) 


Energy  value 


0..10(1) 

0..10(1) 

0..10(1) 

0..90(5) 

0..90(5) 

Calculated 

15. .65(5) 

0..90(5) 

0..90(5) 

15 

0..10(1) 

0..90(5) 

0..90(5) 

15. .65(5) 

0..90(5) 

0..90(5) 

In  progress 

20. .65(5) 

0..10(1) 

0..90(5) 

0..90(5) 

15. .65(5) 

0..90(5) 

0..90(5) 

15. .65(5) 

0..10(1) 

0..10(1) 

0..90(5) 

0..90(5) 

Calculated 

15. .65(5) 

0..90(5) 

0..90(5) 

15 

0..10(1) 

0..90(5) 

0..90(5) 

15. .65(5) 

0..90(5) 

0..90(5) 

In  progress 

20. .65(5) 

0..10(1) 

0..90(5) 

0..90(5) 

15. .65(5) 

0..90(5) 

0..90(5) 

5.1.  Relative  GB  energy  for  specified  GBM  averaged  on  possible  GBIs 

Grain  boundary  energies  were  systematically  calculated  for  various  sets  of  lattice  rotations  (exact 
misorientations)  corresponding  to  given  scalar  misorientations.  Database  for  all  Bunge  angle  combinations 
for  scalar  GBMs  from  1  to  64  degrees  was  created.  For  each  integer  value  of  scalar  misorientation  all 
combinations  that  gives  8  value  in  ±0,02  degree  interval  (0,0,0)  ( (p ,  % ,  6 )  were  derived.  It  was  found  that, 

for  each  scalar  GBM,  arrange  of  three  corresponding  Bunge  angles  consists  from  250  to  830  variants. 
Examples  of  such  arranges  presented  as  tables  5.2  and  5.3. 
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Table  5.2  (8=3  deg) 


N 

9 

X 

0 

1 

0 

3 

0 

2 

1 

1 

3 

3 

1 

4 

0 

4 

2 

2 

3 

5 

2 

5 

0 

6 

3 

3 

3 

249 

90 

87 

0 

250 

90 

90 

3 

Table  5.3  (£  =  29  deg) 


N 

(P 

X 

0 

1 

0 

0 

29 

2 

0 

19 

22 

3 

0 

22 

19 

4 

0 

29 

0 

5 

1 

1 

29 

6 

1 

20 

22 

7 

1 

23 

19 

829 

90 

71 

22 

830 

90 

90 

29 

Then,  for  all  Bunge  GBMs  from  particular  arrange,  GB  energies  were  estimated,  using  software 
developed  in  the  frames  of  the  project.  Two  variants  of  GBI  were  used  with  above  misorientation  arranges, 
namely  (20,  20)  and  (30,  30).  Some  results  for  (20,  20)  GBI  are  presented  on  Figure  5.1. 
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Fig.  5.1.  GB  energies  for  misorientations  of:  (a)  2  deg.;  (b)  3  deg.;  (c)  10  deg.;  (d)  13  deg.;  (e)  29  deg  and 
(f)  45  deg  presented  as  dependence  on  Bunge  angles  variant  number  in  corresponding  dataset 
(proportional  to  first  rotation  angle,  for  GBI  (20,  20) 
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Misorientation,  s 


(a)  (b) 

Fig.  5.2.  Averaged  for  all  orientation  variants  GB  energy  as  function  of  misorientation  (a)  GB  Inclination 
=  (20,  20);  (b)  GB  Inclination  =  (30,  30). 


It  is  clearly  seen,  that  for  a  given  particular  scalar  misorientation  number  GB  energy  may  differ  up 
to  3  times  (Fig.  5.1).  The  calculated  for  GBI  (20,20)  and  (30,30)  values  of  GB  energies  were  averaged  for 
all  misorientations.  It  is  interesting  to  note  that  the  averaged  energy  dependence  on  misorientation  looks 
similar  to  that  usually  used  effective  GB  energies/mobilities  [10]  Figure  5.2. 

5.2.  Database  validation  on  available  experimental  data 

For  specific  GB  geometry  (non-symmetric  tilt  GB  equivalent  to  recent  experimental  study  [5] 
(Fig.  5.3))  GB  energy  was  estimated  as  function  of  GBI,  to  compare  MD  simulation  results  with  available 
experimental  data  (Fig.  5.4). 

It  can  be  concluded,  that  averaged  experimental  and  MD  simulated  dependences  are  not  in  contradiction. 


Grain  B  [001]A  G.B.  Grain  A 


Fig.  5.3.  Bicrystal  with  a  [1  1  0]  asymmetric  tilt  boundary  [15].  The  misorientation  is  defined  as  the  angle 
0  (=  38,9  deg)  between  the  [0  0  1]  directions  of  two  single-crystal  grains,  and  the  inclination  is 
defined  as  the  angle  cp  between  the  boundary  plane  and  the  bisector  of  the  misorientation 
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(a)  (b) 

Fig.  5.4.  Dependence  of  the  tilt  grain  boundary  energy  on  inclination  cp:  (a)  MD  estimated,  lines 
correspond  to  deep  spline  averaging  (dash  green)  and  bin  (20)  averaging  (solid  blue)  and,  (b) 
experimentally  measured  in  [15]. 


5.3.  Comparison  of  different  potential  usage  effect  on  relative  GB  energy. 

Using  advanced  MD  modeling  software  GB  energies  as  functions  of  GBM  and  GBI  was  being 
continued  for  bee  bicrystal  or  general  case  of  not  special  GB  geometry  were  compared  for  cases,  which 
utilizes  essentially  different  interatomic  potentials.  Calculations  done  with  simple  Bom-May er  potential 
were  compared  with  MD  modeling  mns  with  Lennard-Jones  potential: 

<P>  =  ^jCcr/r,)12  -(<r/rs)6},  ^ 

j*i 

(a  =  0.22  nm,  s  =  0.08  eV,  i.e.  a  =  0.28665  nm).  GBI  was  defined  with  GB  normal’s  angles  co  and  §  (see 
Fig.3.18). 

At  first,  GB  energy  was  calculated  for  (0,  0,  0)  (20,  20,  20)  GBM  and  GBI  parallel  to  100  plane  of 
part  of  bicrystal  with  orientation  (0,0,0). 

Using  above  misorientation,  map  for  relative  grain  boundary  energy  as  function  of  GBI,  (normal  to 
GBI  direction  defined  with  (j)  and  oo  angles)  using  angular  step  of  2  deg  for  both  (|>  and  oo  (about  2000 
points)  was  built  and  compared  with  the  same  map  for  Born-Mayer  potential  (Fig.  5.5).  Also  for 
comparison  of  actual  data  GB  energy  section  for  one  particular  4>  (30  deg)  angle  was  plotted  (Figure  5.6). 
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0  20  40  60  80 


Inclination  angle  0 

(c) 

Fig.  5.5.  Relative  GB  energy  for  GBI  (0,  0,  0)  (20,20,20).  (a)  Bom-Mayer  map  (b)  -  Lennard-Jones  map 
and  (c)  Lennard-Jones  map  rescaled  for  4  levels. 


Fig.  5.6.  Relative  GB  energy  for  GBI  (0,  0,  0)  (20,20,20)  calculated  with  Lennard-Jones  (solid  line)  and 
Bom-Mayer  (dashed  line)  potentials. 


STCU 


PROJECT  P-262  (EOARD  068002)-  FINAL  PROJECT  REPORT 


FPR  PAGE  59 


It  is  well  seen  from  the  Fig.  5.5  maps  and  Fig.  5.6  that  peaks  for  GB  energy  calculated  with  different 
potentials  are  very  close  in  positions,  and  these  positions  correspond  to  low-indexed  directions.  It  is 
obvious  that  usage  of  the  Lennard-Jones  potential  results  in  much  sharper  GB  energy  maximums,  however 
there  are  no  essential  differences  in  their  angular  positions  as  compared  to  Bom-Mayer  potential. 

Hence,  it  can  be  concluded,  that  GB  energy  has  its  maximum  values  when  GBI  normals  are  within 
the  certain  crystallographic  zone  of  one  of  the  bicrystals  mainly  due  to  GB  geometry  properties,  and 
potential  selection  affects  the  shape  and  sharpness  of  the  extremums. 

For  the  MC  Potts  modeling  of  GB  migration  implementation  of  any  “sharp”  dependence  of  effective 
GB  mobility,  looks  unreasonable.  Difference  in  effective  GB  mobilities  which  exceeds  1:100  results  in 
maximum  velocity  for  low  mobile  GB  less  than  0,005  MU  per  MCS.  Typical  times  for  grain  growth 
modeling  are  500-1000  MCS,  for  recrystallization  modeling  80-250MCS.  Hence,  just  fixing  of  low- 
mobile  GB  will  occur  within  the  modeling  volume.  Therefore,  usage  of  Bom-Mayer  approach  looks 
preferable.  It  is  also  evident  that  calculations  with  Born-Mayer  potential  provide  more  detailed 
information  on  GB  energy  dependence  on  GBI  and  GBM. 
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6.  p-Ti  RECRYSTALLIZATION  EXPERIMENT  AND  ADVANCED  MC 
MODELING. 

Commercial  P-titanium  alloy  having  actual  composition  Ti  -  1.5  (wt.%)  A1  -  6.8  %Mo  -  4.5%Fe  - 
0.2%O,  (TIMETAL-LCB)  was  selected  for  analysis  of  recrystallization. 

The  initial  microstructure,  texture  (ODF),  spatial  distribution  of  the  stored  energy  of  deformation  are 
accounted  when  modeling  to  compare  how  differences  in  nucleation  mechanisms  and  kinetics  affect 
recrystallization. 


6.1.  Recrystallization  in  cold  drawn  LCB  alloy 


Experimental.  The  thermomechanical  treatment  route  of  LCB  rods  included  solid-solutioning  at 
temperature  of  single-phase  P-field  at  850°C  during  1  hour  with  subsequent  water  quenching.  As  a  result 
single-phase  metastable  P-condition  with  average  P-grain  size  of  about  180  pm  was  obtained. 

Quenched  specimens  were  cold  drawn  with  70%  reduction  to  the  final  diameter  8  mm,  then  cut 
into  specimens  of  60  mm  for  subsequent  bulk  (RHT)  performed  with  direct  electric  current  resistant 
heating,  20  C/s  to  different  peak  temperatures  within  the  recrystallization  interval  (Figure  6.1). 

Two  comparative  Cases  A  and  B  were  run  to  estimate  how  accounting  of  exact  GB  geometry 
affects  recrystallization  in  microstructures  similar  to  those  observed  experimentally  in  LCB  titanium  alloy. 


Modeling.  Modeling  volume  was  constructed  to  mimic  cold  drawn  condition  including  texture 
and  microstructure.  In  both  Cases  3  and  4  nuclei  were  allowed  to  form  only  in  misorientations,  selected 
from  Table  6.1.  Proper  selection  of  nuclei  orientations  allowed  recrystallization  front  motion  be  traced  by 
using  exact  GBM+GBI  approach  (Case  B)  while  in  Case  A  effective  grain  boundary  mobility  followed 
commonly  used  dependence  on  scalar  misorientation.  All  other  initial  parameters  and  conditions  were  the 
same  in  both  cases  A  and  B. 


As  expected,  predicted  recrystallization  kinetics  are  different  due  to  the  slower  migration  of 
unfavorably  inclined  parts  of  RX  front  in  Case  B  (Figure  6.2).  Some  difference  in  recrystallized 
microstructures  are  seen  from  Fig  6.3.  For  Case  4  at  50MCS  (80%  Rex.)  it  is  clearly  seen  preferred  GB 
inclinations  and  grains  arranged  in  vertical  lines  (along  the  long  axes  of  the  deformed  grains). 

Similar  regular  grain  arrangements  are  seen  on  the  experimental  LCB  grain  boundary  maps  for  95% 
recry stalized  material  (Fig  6.4). 


Fig.  6.1.  Optical  (left)  and  EBSD  (right)  images  of  a)  solid  solutioned;  b)  cold  drown;  c)  50% 
Recrystalized  and  d)  100%  recrystalized  material. 
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Fig.  6.2.  Predicted  recrystallization  kinetics  for  Cases  A  and  B. 


50MCS 
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100  MCS 


Fig.  6.3.  Simulated  microstructures  for  Cases  1  (left)  and  2  (right).  GB  color  is  proportional  (red 
GBs  >15  deg  mis,  other  <15,  blue  <  2  deg) 


Fig.  6.4.  LCB  95%  Rex.  EBSD  Grain  Boundary  Map  (400*500  pm)  shows  grain  arrangements  (marked 
with  green  color)  which  looks  similar  to  Case  4  simulated  patterns  (right). 


6.2.  Recrystallization  in  specially  textured  3D  MC  volume 


For  15  selected  grain  orientations  GB  energies  were  calculated  as  those  depending  on  GBI:  full 
maps  of  GBI  dependant  GB  energy  were  built  for  grain  misorientations  listed  below  (Table  6.1)  with  one 
grain  in  (0,0,0)  position,  as  well  as  for  possible  pair  misorientations  from  the  list.  Calculated  maps  are 
shown  on  Figs.  6.6  and  6.7.  Using  these  orientations,  an  artificial  as-deformed  microstructure  was 
constructed,  which  was  used  for  MC  simulation  with  taking  into  account  exact  grain  boundary  geometry 
(GBM+GBI).  Such  method  of  construction  for  initial  as-deformed  3D  modeling  volume  guarantees 
presence  of  the  finite  set  of  possible  misorientations  in  modeling  volume,  with  known  GBI  dependencies 
for  GB  energy.  Effective  GB  mobilities  were  assumed  to  be  proportional  to  corresponding  GB  energy. 
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Table  6.1. 


N 

cp ,  deg 

x>  des 

6 ,  deg 

1 

21 

21 

3 

2 

30 

30 

3 

3 

0 

0 

3 

4 

3 

2 

15 

5 

0 

0 

15 

6 

31 

22 

12 

7 

32 

44 

9 

8 

45 

45 

0 

9 

10 

35 

0 

10 

12 

37 

0 

11 

14 

39 

0 

12 

0 

0 

30 

13 

18 

48 

0 

14 

37 

7 

0 

15 

45 

45 

45 

so 


20  40  €0 

Inclination  angle 


(a)  (b) 

Fig.  6.  6.  Dependence  of  relative  GB  energy  on  GBI  for  high  symmetric  GBMs:  (a)  (45,45,45)  and 
(b)  (45,45,0). 

Two  MC  Recrystallization  simulation  runs  were  completed,  with  (case  2)  and  without  (Case  1) 
accounting  energy  GBI  dependencies.  Starting  from  the  same  as-deformed  microstructure  presented  on 
Fig.  6.5  recrystallization  starts  with  300  nuclei.  Orientations  of  nuclei  were  selected  randomly  from  Table 
6.1.  As  expected,  due  to  the  reduced  average  mobility  when  GBI  dependencies  are  accounted,  Case  2 
performs  slower  kinetics  of  recrystallization  (Figure  6.8).  Microstructure  evolution  for  Cases  1  and  2  is 
presented  on  Fig.  6.9. 
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Fig.  6.7.  Dependence  of  relative  GB  energy  on  GBI  for  general  GBMs  from  table  6.1. 
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Fig.  6.8.  Recrystallization  kinetics  for  Cases  1  and  2. 


30MCS 


30MCS 


40  MCS 


40  MCS 
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Fig.  6.9.  Microstructural  evolution:  (a)  Case  1;  and  (b)  Case  2. 


Microstructures  at  60  MCS  for  Cases  1  and  2  show  some  difference  in  GBI  population,  which 
random  for  Case  1  and  partially  follows  some  definite  planes  for  Case  2.  (Figure  1 .39). 


Fig.  6.10.  Some  favorite  GB  planes  for  Case  2,  60  MCS. 
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New  program  for  MC  Potts  modeling  of  microstructural  and  textural  evolution  within  the  3D 
modeling  volume  was  tested  and  installation  package  on  CD  was  generated.  Project  software  CD  was  send 
to  partners  (Dr.  Wynn  Sanders,  EOARD  and  Dr.  Lee  Semiatin,  WPAFB). 

Two  screenshots  are  presented  on  Fig  6.1 1.  Installation  package  is  included  with  project  CD. 


(a) 


(b) 


Fig.  6.11.  MC  Potts  Rex.  and  GG  software:  (a)  modeling  volume  analysis;  (b)  searching  for  GB  energy 
value  within  included  database. 


Distinguishing  features  of  the  software  are: 

Exact  GB  description  into  the  3D  MC  code  was  implemented.  Now  usage  of  an  exact  GB 
misorientation  (five  fold  GBM+GBI),  instead  of  the  usual  scalar  approximation,  is  possible; 
Procedure  for  local  GBI  estimation  according  to  developed  procedures  (see  section  1); 

Software  allows  including  of  the  GB  properties  database,  generated  by  GB  MD  simulation  routine 
for  instant  usage  during  MC  Potts  modeling  run; 

Manual  setup  for  grain  orientations  in  the  “as  deformed”  state  now  is  possible. 
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7.  RELATIVE  GB  EFFECTIVE  MOBILITY. 

To  estimate  GB  effective  mobility  for  3D  MC  simulation  technique,  it  is  also  needed  to  determine  the 
mobility  of  a  single  GB  as  function  of  both  GBI  and  GBM,  and  then  calculate  relative  GB  mobilities  as 
functions  of  GBI  for  each  particular  GBM. 

Two  GB  geometries  were  selected  for  MD  simulation  of  GB  mobility  as  function  of  GMB+GBI.  The 
first  one  is  traditional  bicrystal  half-loop  geometry.  This  geometry  allows  to  investigate  the  shape  of  a 
moving  grain  boundary  and  use  it  as  a  source  of  data  with  regard  to  the  GBI  dependence  of  GB  mobility 
both  experimentally  and  with  MD  simulations. 

For  given  misorientation  (Figure  7.1a)  half-loop  geometry  allows  to  estimate  GB  effective  mobility  for 
all  possible  GB  inclinations.  When  GB  migrates,  its  shape  is  changing  towards  the  polygon,  with  sides 
close  to  the  inclination  of  the  slowest  GB  segments.  In  modeling  experiment  with  GBI  dependent  GB 
velocities  of  local  GB  segments  will  differ  from  reference  Vi  and  V2  velocities  for  GBI  1  and  GBI  2, 
correspondingly  (Fig.  7.1a). 

By  using  quite  similar  (quarter-loop,  Figure  7.1b)  geometry  the  above  was  shown  experimentally  for 
various  metals,  for  instance,  A1  [11]  (Figure  7.2).  On  figure  6a  grain  boundary  segmentation  due  to  the 
differences  in  GBI  dependent  mobility  is  clearly  seen. 

However,  such  approach  requires  a  huge  amount  of  computational  time.  As  a  result,  in  [3]  simplified 
“molecular  statistic’  approach  was  used  instead  of  MD,  and  just  for  2D  case.  Hence,  we  plan  to  study  only 
several  key  misorientations  with  this  approach. 


(a)  (b) 

Fig  7.1.  Steady  state  GB  migration  (half-loop)  geometry  as  proposed  for  MD  estimation  of  GB  mobility  (a) 
and  (b)  -  geometry,  often  used  for  experimental  GB  mobility  estimation. 
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Fig.  7.2.  Experimental  GB  mobility  estimation  for  A1  alloy  [1 1].  (a)  low  temperature,  essential  GBI 
dependence  (three  nearly  flat  segments),  (b)  high  temperature,  almost  no  GBI  dependence. 


Generally,  the  driving  force  for  grain  growth  is  the  reduction  of  the  free  energy;  this  reduction  is 
accomplished  by  the  migration  of  the  grain  boundary.  For  half-loop  geometry  the  source  for  the  driving 
force  is  the  area  of  the  boundary,  however  other  sources  can  contribute  to  the  driving  force  and  hence  to 
grain  boundary  motion.  Such  sources  can  be  given  by  a  gradient  of  any  intensive  thermodynamic  variable, 
for  example,  a  gradient  of  temperature,  pressure,  a  magnetic  field  (Figure  7.3)  [12-14],  etc. 

Hence,  by  introducing  the  external  driving  force  of  any  nature,  not  depending  on  GB  curvature,  it  is 
possible  to  study  GB  effective  mobility  for  a  simple  flat  GB,  having  the  same  GBI  in  each  point.  For 
modeling  purposes  no  matter  what  is  the  driving  force  nature.  It  just  can  be  introduced  as  constant  vector 
adding  to  the  interatomic  interaction  potential. 

Initial  calculations  were  started  to  show,  that  energy  effect  both  from  additional  driving  force  and  GBI 
can  characterize  the  GBI  effect  on  GB  effective  mobility. 


Fig.  7.3.  Example  of  spatial  orientation  of  the  rolled  sample  with  respect  to  the  magnetic  field  H,  provides 
essential  driving  for  GBs  force  due  to  the  anisotropy  of  the  magnetic  susceptibility  of  titanium 
(magnetic  annealing  of  titanium,  [12]) 

Two  techniques  were  initially  developed  using  previously  proposed  approaches  for  estimating  the 
mobility  of  a  single  GB  as  function  of  both  GBI  and  GBM. 
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The  first  one  is  an  attempt  to  model  GB  migration  in  traditional  bicrystal  half-loop  geometry 
(Fig.3.1a)  This  geometry  allows  investigating  the  shape  of  a  moving  grain  boundary  and  using  it  as  a 
source  of  data  with  regard  to  the  GBI  dependence  of  GB  mobility. 

The  second  technique  for  MD  GB  mobility  estimation  is  based  on  the  introducing  of  external  driving 
force  which  does  not  depend  on  GB  curvature.  This  technique  was  being  developed  to  the  study  of  GB 
effective  mobility  for  a  simple  flat  GB  having  the  same  GBI  in  each  point. 

First  set  of  calculations  was  completed  for  the  (0,0,0)  (20,20,20)  misoriented  boundary.  The  Born- 
Mayer  potential  was  used  with  adding  of  as  constant  vector  force  for  affecting  each  atom  within  the  thin 
volume  close  to  GB  (within  3  lattice  periods  from  GB). 

For  this  particular  case  no  essential  influence  an  external  force  on  GBI  energy  dependence  was  seen. 
The  only  effect  was  nearly  constant  difference  in  absolute  values.  The  normalized  dependencies  were 
similar  (see  Fig.  7.4).  It  is  also  possible,  that  we  selected  GB  area  too  wide,  and  as  the  result  could  not  see 
differences  in  the  “relaxed”  energies.  So,  we  apply  additional  driving  force  for  thinner  layer  around  GB. 


Fig.  7.4.  Relative  GB  energy  for  (0,0,0)  (20,20,20)  misoriented  GB  as  function  of  GBI  component  c|).  Open 
circles  are  GB  energy  values  when  additional  driving  force  is  effective. 


Fig.  7.5.  Example  of  flat  relaxed  GB  in  MD  simulation 
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An  attempt  to  model  GB  migration  in  traditional  bicrystal  half-loop  geometry  was  performed.  This 
geometry  allows  investigating  the  shape  of  a  moving  grain  boundary  and  using  it  as  a  source  of  data  with 
regard  to  the  GBI  dependence  of  GB  mobility. 

First  two  runs  for  Z9  GB  and  5  runs  for  general  symmetric  tilt  boundaries  were  completed.  It  was  shown, 
that  the  movement  of  such  GBs  in  modeling  volume  are  too  slow,  and,  as  expected,  they  require  a  huge 
amount  of  computational  time  to  build  even  simple  one-dimensional  graph. 

Hence,  to  show  effects  of  GBI  inclination  accounting  on  grains  shape  an  recrystallization  kinetics,  it  was 
decided  to  use  data  for  relative  GB  energy  as  effective  mobilities  for  3D  MC  Potts  modeling  code. . . 


Fig.  7.6.  Migrating  hemispheric  GB  geometry  used  in  MD  simulation. 
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8.  TOOL  FOR  ANALYSIS  OF  NOT  RELAXED  GB  GEOMETRY 


As  additional  tool  was  developed  (Fig.  9.1)  for  analysis  of  not  relaxed  GB  geometry.  It  allows  setting 
lattice  misorientation  in  terms  of  rotation  axis  and  rotation  angle  and  visualizing  the  CLS  pattern  by 
showing  sites  of  the  misoriented  lattice  which  occupy  the  same  or  slightly  different  (difference  in  position 
not  more  than  chosen)  positions. 


(a) 


*  *  ^  ^  £  £  XX 

X^x.  X  «s>  ^  ^  ^  ^  ^  ^ 

O WX'V*  \w!b  X  X  <*  *4  *4  *4  ^  ^  ^ 

Vi'Sk <» <* <» 5* 54 Cj ^ ^ 5i ^ XXvXXXX 
^J^?*?****!*  5*  5*  5*^^^  ViV^V 

fi  fi  55  5*  5*  4  4  4 

^V^ivlMvs^S^fssi  fj  fj  4  5$  5$  55  % 

ji  \js *s ^ ^ ^  **  f*  &  $> 

'  *  *  *VMNM*jM*M**ii  X$S$xSxSxSi*«^i^.^.^  ^  § . 


4  §  0  §  i  §  SHMMWMW*»  \*SSx $ . 

l>  s>  3  <5  0  8  8  ***<M*!*sk>  Vsss*s*s*s^^^^^ . 

X*S*s$s$s$*§*§*§i£i^ 

S?^<V*V^  X  $  I*  ^  ^  ■&****<&*  \\*s*x$x$x$xSxSxS\l\ 

<i  <1  |5  ?5  $  $  $  \vfS*s$x$x§xlx$x$\ 

f  *S  \  5  \  5  \  ^  ^  ^  ^  ^  ^5 

3^  3)  ?*  ?  ?* 

'  . ~  S><i 

>4  Xc> 


(b) 


.  'ktWiA.Vi  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  OOO*e0OOOOOOOOOOOOOOOOOOOO 

O  °  0  0)^^00^00  oooooooooooooo 
ri  °  °  OQ30DioO  T)  (ti  CD  'o  d  (5  (5  o  oooooooooo 

Ta  ^  ^  O  O  OPCD^O°C^^DOCo^OCbcticO0Di^J(tiCD(Sl(f'  o  O  O  O  O  O 

\ ■SSa^rf^ $$ Si^r^ $$  ^S^5*iXiri<rs. %  O  O  ®  O  Q  ^CP  03  ePeP  CO  (^C^cfloGDc^OTD  0D  00^0,0  O 

^^5*55*5^^  ooO(jgoSc0g(?(Po0  g  cP  cP  cto  cr  cPcdcdc^o0^  o 

\^SsSs:i^i^  ^^^^SSsSsti^i^i^^SKSSs  O  OfttP^S  0  5(S(J  8  9  O  Q  Q  8  S  0  0  (?  <?  <P  O  Q  tf  (P  O 

SsSi:i^i|,i  55  o  cAcocWpftdOcPfpS.  9)  05(55,0  BfiOijjS  0  0  (3  fl  ff  O  , 

AXX*X  5i ** X  O  ipixiiH-, '^'o'-'ini' cn q-, o S:. '+0:^ :, ' o % Qjcorfj S  Pi  drD.v'S  ?.  o  ,: 

55  ?A  Oeoeatffi)C%Q^  ^cpcocb^o^tccyo^ PjO  q 

>^.'i ^ ^^5 A>r^  x ^ ^ oQoaj8eoeJeppef^(D(&§)J(p(D(i(t^cP'o 

oPjcutp  jd  %  9)  a)  <5  J3  6O3go8O0g(f0o0Q(9(£)o 
cflcPcociro^^cocpJioPjQjcj)^?)  P)  OagiSQog  o  S  9  o 

,a  „  _  ,  _ T.  .._  -  ^  cPcPQja  ^cPcocfedo0oPooci-p'°octd3Cj]Ho9)Cl)a)(prp  ft  9)  o 

^  ^5  ^  \^S:ES^^  %  rV^f  ^ =S! Si^r^  0Co  g  ocPueo^cPcDC&dP  cPcPdjd,  d,  oPqPqoc^o^  Pjo 

^  X-X 1  ^  ^  i  ^  ^  r^.  s  1 5  ^  0dcS(?S  0  0|JD  5^00  0  ^PCg^CpCpQj^  cPcPoO 

ri  X^i^i':ir;^^^|5g'^sss:  Pjcscjj^jS  PxScP^BQctiSo  0  QogSfiPo  g  J(?(D 

ri  ^  °CpQ0Cl30°0CtPQoQw0'%Q>Q0cprp  ^^CE^SSdcSflS  9  0 

^  ^  °,V  ,D  c-,  |£ i+'iS'  □:,  ■-;-,  cPrPtf1  Co O  O Cb 'dccto.-.o  %  CnrDgj  Jli  ft  d 

^  X  -  \x  xK^W  °0oeP^0F®®  .?Xd'  CD  d,  a°o°cPcoq5°oft^ 

XXiX-'XX  ^  °1o)0)£i09ogf0O0gff0o0Q  cP<9  (D  dPcPcPo 

;§XX"X% 5i ^  VSXXX>  :  '>.co1Sixft,d,a,1?;5 ftd^SesSSsQoeffleo 

/!- s|!s  b^lt^W£^ 1 

X  X  x  ^  ^  ^  "x  ^  X  ^  ^  XX  X  V  °  (ScPt?  5  0O(j(j8  0O(5DSflc0ji'(9i)(j)0  cPcPo  o 

\  !51'“  — \  00  0  »«-0o0„"/".0.0„"»-°o0.0„»?o0 

(d)  (e) 


000*0000000000000000000000 
o  o  oa>tf>  £  9  05.00000000000000000 
O  O  COCbC0O<boD0D  ft  ftcp^l  0000000000000 
O  O  <£>150  Cp<D<fcO  ePooegO  ftODtfJ  JD  OOOOOOOOO 


ocP<5^cPcocb(XoCboo5o95CboD0DftQi<5)<f  9ogo0oc 
0  0  0  g  Pc  ©  gcPcDi^^cPQjc^o^joc^D^ftcrcpS^ai  c 
®  ^0555  0  0(5^900  01900  (PiP^iXcPo^cftc^oPooo  o 
OooqjO ftODco^  ftd,p  S  9  o  51  8  9  0  0  8  &  ©  Q  cPrno  o 

0  «D  (5  ^*cPojC^OOCOCt>C0OCboDcp  %  ^.(D(P  000^  9  0  0  0 

°  o  q  8  6  ©  0  cPp®  i^cPcP^o  o®ooc0oftcbcD(^>ftdi  o  o 
°  a>(j)  o9.aig8  9  <5  g  8  Q  ©  §8  id  <5  tf>cP<Dq>o°cao  o  o 
5  ochc^o9oocx-jo'55Q)cP(-p  0O(jiD  9  Ofl  oSofl  cf  5*o  o  o 
°  O0T0  O  n  fo  on  5o  o  „  fto  o  °  0  o 


Oo°«o 


°o°°o 


?ooo°o°;#c 


.0°on0ooo r  °o 


(6 


Fig.8.1  CLS  visualization  tool  (a)  and  CLS  patterns  changing  from  LABs  to  HABs  (twist  GB,  0=  2  deg. 
(b),  4  deg.  (c),  6  deg.  (e)  12  deg  and  (f)15  deg). 
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2.  Shevchenko  S.V.,  Ivasishin  O.M.,  Gazder  A.,  Pereloma  E.,  3D  Monte-Carlo  simulation  of 
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1)  MC  -  Monte  Carlo 

2)  MU  -  Modeling  unit,  modelling  units  form  the  Monte  Carlo  Potts 
modelling  volume . 

3)  MD  -  molecular  dynamics ; 

4)  GB  -  grain  boundary ; 

5)  GBI  -  grain  boundary  inclination ; 

6)  GBM  -  grain  boundary  mi s orientation ; 

7)  M0(x,  y,  z)  -  Current  MU  position . 

8)  R  -  GB  curvature  radius . 

9)  n  -  normal  for  planar  linear  approximation 

10)  m  -  GB  tangential  plane  normal 

11)  M  —  GB  misorientation  (matrix  form) 

12)  8  -  GB  misorientation  (scalar  form) 

13)  ODF  -  orientation  distribution  function 

14)  TTC  -  tilt/twist  component 

15)  BW  -  texture  band  width; 

16)  SBG  -  subgrain  size ; 

17)  LAB  -  low  angle  GB; 

18)  HAB  -  high  angle  GB; 

19)  MCS  -  Monte  Carlo  step; 

20)  CA  -  Cellular  automata ; 

21)  STGB  -  Symmetrical  tilt  grain  boundary; 

22)  VMD  -  Visual  molecular  dynamics ; 

23)  CLS  -  Coinside  Lattice  Structure  (Sites) 
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